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systems. 
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Rehevot,  Israel. 
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DESIGN  OF  A  MULTIPLE  INPUT-MULTIPLE  OUTPUT  FLIGHT  CONTROL  SYSTEM 


CONTAINING  UNCERTAIN  PARAMETERS 

I.  INTRODUCTION 

The  approach  used  to  design  the  flight  controller  in  this  research 
project  is  the  Quantitative  Feedback  Theory  (QFT)  technique  developed  by 
Dr.  Isaac  Horowitz  of  the  Weizmann  Institute  of  Science,  Rehovot, 

Isreal.  He  used  this  technique  on  several  projects  for  the  U.  S.  Air 
Force  (1A).  This  frequency  response  technique  was  used  to  successfully 
design  controllers  for  several  difficult  systems;  for  example,  the 
FY-16  CCV  aircraft  (12). 

The  U.  S.  Air  Force  is  interested  in  this  technique  as  a  method  to 
synthesize  multiple  input -multiple  output  (M1M0)  controllers.  In  this 
thesis,  a  lateral  flight  controller  is  designed  for  the  Flight  Propulsion 
Control  Coupled  (FPCC)  aircraft. 

1.  Background 

An  aircraft  sitting  on  the  ground  is,  in  general,  a  stable,  decoupled 
system.  In  the  air,  an  aircraft  becomes  a  nonlinear  and  highly  coupled 
M1M0  system  whose  model  contains  parameter  uncertainty.  The  aircraft 
system  Includes  many  inputs,  some  from  flight  control  surfaces,  e.g. 
rudder,  ailerons,  spoilers,  flaps,  elevator  and  canards.  Also,  inputs 
come  from  propulsion  units  or  engines  (e.g.  thrust),  and  from  other 
sources  such  as  gusts.  Some  of  the  system  outputs  controlled  are 
velocity,  angular  displacements  and  angular  rates  of  the  aircrait  body 
relative  to  some  reference  axis  system. 


Ideally,  each  Input  affects  only  the  one  or  two  outputs  it  is 
designed  to  affect.  In  practice,  one  input  affects  many  outputs  because 
the  inputs  and  outputs  of  an  aircraft  system  are  strongly  interrelated, 
or  coupled.  The  input  to  output  relationships  usually  are  not  linear 
nor  are  they  known  with  much  certainty.  For  example,  the  ailerons 
control  the  roll  rate,  which  controls  the  bank  angle.  But,  an  aileron 
deflection  also  causes  a  yaw  rate  and  a  pitch  angle  change.  The  deflected 
surfaces  cause  asymmetric  drag  on  the  wing,  thus  a  yaw,  and  the  changing 
lift  vector  resulting  from  the  change  in  bank  angle  causes  a  pitch 
change. 

The  classical  mathematical  model  considers  the  M1M0  relationships 
invariant  and  linear  in  a  small  region  around  an  equilibrium  point 
(6:152-155).  Next,  the  MIMO  system  is  decoupled  into  several  independent 
single  input-single  output  (SISO)  problems  (5:162).  The  designer  learns 
the  art  of  decoupling  mostly  through  experience.  Only  experience  helps 
decide  which  simplifications  are  appropriate.  The  resulting  designs  are 
often  successful,  but  only  after  much  trial  and  error  (14:1). 

Currently  many  controller  designs  use  modern  or  optimal  control 
techniques.  These  control  design  techniques  use  the  state  space  repre¬ 
sentation  of  the  system  model  (5:523-601).  These  techniques  normally 
require  many  complex  matrix  operations.  Robustness  using  these  design 
techniques  is  not  guaranteed  and  is  difficult  to  obtain  when  unavailable 
state  information  must  be  estimated.  (See  Reference  11  for  more 
critique).  Robustness  is  defined  in  this  paper  as  a  single  controller 
which  maintains  the  desired  aircraft  performance  throughout  the  flight 
envelope.  Another  method  studied  at  the  Air  Force  Institute  of 


Technology  (AF1T)  and  In  the  Flight  Dynamics  Laboratory  (FDL)  at 
Wright-Patterson  Air  Force  Base,  is  the  output  feedback,  high  gain, 
asymptotic  method  of  Professor  Brian  Porter  of  the  University  of 
Salford,  U.K.  This  method  guarantees  robustness  in  the  design,  but  the 
amount  of  robustness,  or  range  of  uncertainty  over  which  the  design 
performs  satisfactory,  is  not  predictable. 

Another  method  developed  by  Dr.  Horowitz  (10-17)  and  used  in  this 
thesis,  guarantees  to  satisfy  the  system  performance  specifications. 
System  responses  to  inputs  or  plants  which  are  within  the  expected  input 
or  plant  ranges  are  guaranteed  acceptable.  Thus,  the  design  satisfies 
the  system  specifications  of  many  plant  and  input  combinations  by 
designing  for  only  a  finite  number  of  combinations.  Therefore  the  QFT 
technique  is  robust.  Dr.  Horowitz's  technique  begins  the  design 
assuming  uncertainty  in  the  plant,  the  disturbances  and  the  performance 
specifications.  The  differential  equations  or  the  state  space  model  for 
the  system  represents  the  MIMO  plant.  Several  equivalent  S1S0  systems, 
which  contain  disturbance  and  parameter  uncertainty,  are  synthesized 
through  elementary  matrix  operations.  This  reduction  from  a  MIMO  system 
to  several  SISO  systems  is  mathematically  rigorous.  A  design  is 
performed  for  each  S1S0  system  using  the  system  performance 
specifications  and  information  derived  from  the  MIMO  system.  If  the 
designed  controller  synthesized  for  each  SISO  system  meets  all  of  its 
specifications,  then  the  final  controller  is  guaranteed  to  meet  the  MIMO 
specifications  (14). 

The  flight  controller  designed  in  this  thesis  is  for  a  "paper" 
aircraft.  This  hypothetical  aircraft  was  developed  under  a  FDL  contract 
by  Lockheed-Calif ornia ,  Honeywell  and  Pratt-Whitney  (7,8).  This 


hypothetical  aircraft  is  statically  unstable  and  incorporates  canards 
and  jet  flaps  for  maneuverability  (7,8).  The  contractors  provided  FDL 
with  a  computer  simulation  of  this  aircraft,  FPCCS1M  (9).  This  program 
simulates  the  FFCC  aircraft  at  different  flight  conditions.  Additional 
computer  programs  used  during  this  thesis  are  TOTAL  (an  interactive 
design  program),  a  program  to  invert  three-by-three  transfer  function 
matrices  (Appendix  F)  and  programs  generated  by  the  author  of  this 
thesis  specifically  for  the  QFT  technique  (Appendix  F) . 

2.  Problem 

This  thesis  presents  a  design  of  a  MEMO  feedback  lateral  flight 
controller  for  an  aircraft  using  the  Quantitative  Feedback  Theory  (QFT) 
technique  developed  by  Professor  Horowitz.  The  controller  is  designed 
for  the  Flight  Propulsion  Control  Coupled  (FPCC)  "paper"  aircraft,  using 
results  from  FPCCSIM  (a  computer  simulation)  for  the  aircraft  model. 

The  resulting  design  is  compared  qualitatively  to  a  design  previously 
accomplished  for  this  aircraft,  using  Professor  Porter's  state  variable 
MIMO  feedback  technique  (2).  Also,  the  cost  for  achieving  desired 
performance  despite  uncertainties  is  discussed.  The  design  method  is 
outlined  such  that  it  can  be  applied  to  similar  problems  by  anyone  with 
a  background  in  control  theory  and  frequency  response  methods. 

3.  Assumptions 

The  basic  modeling  assumptions  are  a  flat  non-rotating  earth 


reference  frame,  a  quasi-static  airflow,  a  rigid  aircraft  and  a  constant 
aircraft  mass  (4:21-23).  These  assumptions  limit  the  problem  to  one 
achievable  in  one  thesis  and  not  by  the  capability  of  the  QFT  technique. 


the  M1M0  design  approach  in  Chapter  III.  Next,  in  Chapter  IV,  the  FPCC 
aircraft  and  FPCCSIM  computer  simulation  are  described.  The  computer 
program  FPCCSIM  generates  the  aircraft  model.  After  completion  of  this 


this  method  by  the  author  ot  this  paper  (Appendix  F)  and  the  computer 
aided  design  program  TOTAL.  Evaluation  ot  the  design  is  continuous. 
During  the  design  phase,  changes  or  compromises  are  made  in  the  design 
as  appropriate.  This  demonstrates  the  method's  "transparency"  or 
ability  to  evaluate  difficulties  in  the  system  during  the  design  phase. 
The  designer  knows  when  a  problem  exists  and  what  needs  modification  to 
alleviate  the  problem.  The  cost  in  bandwidth  of  any  change  is  evaluated 
quickly.  These  changes  do  not  affect  the  original  assumptions  (see 
Section  1-3)  nor  the  system  performance  specifications  of  the  final 
controller. 

Chapter  Vll  contains  the  computer  simulation  and  compares  the 
results  to  the  system  performance  specifications.  Then  the  design  is 
compared  qualitatively  to  the  previously  designed  controller  using 
Professor  Porter's  method.  A  discussion  and  conclusions  follow  in 
Chapter  VIII.  Recommendations  for  further  studies  using  the  QFT 
technique  are  also  included. 


II.  QUANTITATIVE  FEEDBACK  THEORY: 
SINGLE  INPUT-SINGLE  OUTPUT  (S1SO)  SYSTEM 


1.  Introduction 

This  section  presents  an  overview  of  the  QFT  technique,  as  used  In 
this  thesis.  The  technique  description  is  adapted  from  Chapters  11  and 
111  of  Robert  Betzold's  thesis  (3).  The  explanations  and  examples  are 
essentially  the  same,  with  some  modification. 


2.  General  Procedure 


This  section  presents  the  general  procedure  for  this  technique.  A 
more  detailed  explanation  of  the  technique  follows  in  this  and  the  next 
chapter.  Chapters  V  and  VI  contain  the  design  performed  using  this 
technique  in  detail. 

With  a  single  input-single  output  (SISO)  system  defined,  begin  the 
Quantitative  Feedback  Theory  (QFT)  technique  with  the  specifications  for 
the  system.  The  QFT  technique  requires  the  system  specifications  in  the 
frequency  domain  (Section  11-4).  Next,  the  plant  transfer  function 


matrix  model,  ]?  ■  [p^  ]  ,  is  generated  because  the  QFT  technique 


requires  the  frequency  responses  of  the  plant  (Section  11-5). 

The  system  specifications  and  the  plant  models  help  define  the 
limitations  or  bounds  for  the  system  at  every  frequency  of  interest,  <o, 
(Section  11-8).  Due  to  uncertainties  in  P,  P  has  different  values  at 


each  co.  This  resulting  set  of  P^joo),  forms  a  template  of  p^.  Plant 


templates  are  made  for  a  number  of  frequencies  in  the  region  of  interest 
The  plant  templates  and  the  system  limitations  generate  design  bounds  on 
the  Nichols  Chart  (Section  11-9).  Another  design  bound,  the  Universal 


High  Frequency  (UHF)  Bound  maintains  a  safety  margin  for  stability 
(Section  11-10).  The  Nichols  Chart  [log  (L)  with  loci  of  constant 
magnitude  and  phase  of  L/(l  +  L) ]  is  convenient  for  finding  the  bounds 
on  a  nominal  loop  transmission,  L^.  The  loop  transmissions,  or  open 
loop  transfer  functions,  around  the  M1M0  plant  are  the  basic  means  of 
achieving  the  desired  response  specifications  despite  uncertainties. 
is  designed  from  the  design  bounds  using  the  plant  templates  (Section 
11-11).  From  Lq,  the  compensator  is  synthesized  using  Lq  *  Pq  G  where 
G  is  the  compensator,  is  the  nominal  plant  and  Lq  is  the  nominal  loop 
transmission  (Section  11-12).  With  the  loop  characteristics  known  the 
prefilter,  F,  is  designed  (Section  11-13).  This  completes  the  design 
technique  for  a  S1S0  system. 

3.  Problem  Definition 

Consider  a  single  input-single  output  (S1SO)  system  for  simplicity. 
The  plant  transfer  function,  P,  contains  uncertainties  in  its  parameters 
(gains,  zeros  and  poles).  The  system  specifications  define  the  desired 
response  of  the  system,  to  command  inputs  and/or  disturbances.  There¬ 
fore,  the  problem  is  to  achieve  the  desired  response  despite  the 
uncertainties  in  the  plant. 

In  the  S1S0  system  in  Fig.  2-1,  r(t)  is  the  command  input  and  d(t) 
is  the  disturbance  input,  both  of  which  contain  uncertainties.  P  is  the 
plant,  which  contains  uncertainties.  G,  the  compensator,  and  F,  the 
prefilter  are  designed  to  force  the  system  to  have  acceptable  responses. 
The  system  output,  y(t),  is  a  member  of  the  set  of  acceptable  responses 
Y(t).  The  plant  input  is  x(t). 


Fig.  2-1  Single  Input— Single  Output  (SISO)  Control  Loop  with 
Two  Degrees  of  Freedom 


The  quantities  r(t)  and  y(t)  are  assumed  measurable  and  y(t)  is 
assumed  available  for  feedback.  Since  both  signals,  r(t)  and  y(t)  are 
accessible,  a  two  degree-of-freedom  structure  is  used  (see  Fig.  2-1). 
Thus,  the  designer  has  two  compensator  elements,  F  and  G  (14:13).  Als 
assume  r(t),  y(t)  and  p(t),  where  y(t)  -  p(t) *  x(t)  ,  are  all  Laplace 
Transformable  functions  (14:8). 

Figure  2-1  ha6  four  transfer  functions  of  interest.  The  loop 
transmission,  L,  is  defined  as  L  ■  GP  .  The  system  control  ratio  due 
to  the  command  and  disturbance  inputs  are  respectively; 


TR  -  Y (s) /R(s) 


Td  =  Y(s)/D(s) 


FG  P 
T  +  GP) 


FL 

TTT) 


(  1  +  GP )  (l  +  I.) 


(2-1) 


U-2) 


and  the  transfer  functions  relating  the  plant  input  to  the  command  and 
disturbance  inputs  are  respectively; 


IR  *  X(s)/R(s)  -  +  GP)  (i  +  l)  12-3) 

— G  _T  /  p 

1D  -  X(s)/D(s)  -  Q  +  GP)  “  (1  +  L)  (2-4) 


The  design  specifications  can  impose  constraints  on  any  combination 
of  these  transfer  functions,  but  for  this  example  and  for  simplicity, 
only  the  first  two  functions  are  constrained. 

4.  System  Specifications 

The  system  specifications,  or  closed-loop  response  tolerances, 
define  the  acceptable  frequency  domain  output  response  to  inputs; 
disturbance (s)  in  terms  of  an  upper  limit  and/or  command(s),  in  terms  of 
an  upper  and  lower  limit.  Any  response  within  these  bounds  is  assumed 
acceptable. 

The  response  specifications  are  required  in  the  frequency  domain 

for  this  technique.  Often,  the  system  specifications  are  defined  in  the 

time  domain,  using  characteristics  such  as  maximum  peak  value  (M^) , 

settling  time  (t  ),  time  to  peak  value  (t  ),  final  value  (F  )  and 
s  p  v 

maximum  gain  allowed  (K  ).  These  characteristics  are  based  upon  a 

m 

specific  forcing  function,  such  as  a  step  input.  The  specifications 
might  be  defined  as  a  bounded  region  in  the  time  domain,  where  the 
desired  system  response  to  an  input  is  in  the  region  enclosed  by  the 
upper  (Tjj)  and  lower  (T^)  limits,  as  ir.  Fig.  2-2.  Thu  system  response 
to  a  disturbance  input,  T^,  is  normally  defined  as  | y  ( t ) }  <  |  | . 


Fig.  2-2  Design  Limits  for  a  Step  Input  (Time  Domain) 


To  use  the  QFT  technique,  translate  the  time  domain  specifications 
into  the  frequency  domain.  Do  this  by  generating  a  transfer  function 
which  yields  the  same  response  characteristics  as  the  time  domain  system 
performance  specifications  given.  One  method  is  to  use  pole-zero 
placement  as  described  in  Section  12-2  of  Reference  5.  A  third  order 
control  ratio  model  with  one  zero  is  suggested  for  Lhe  response  to  a 


step  input: 

A(s  +  z. ) 

M  - - - - 

1  T 

(s~  +  2  .i  s  +  •“)  (S  +  r  .  ) 


(2-5) 


The  pole-zero  pattern  of  the  above  transler  function  is  in  Fig.  2-3. 
Adjust  the  locations  of  the  roots  until  the  response  of  the  modeled 
control  ratio  matches  the  desired  upper  and  lower  bound. 


.m2 


ii 


During  the  modeling  of  the  equivalent  frequency  domain 
specifications,  consider  the  frequency  domain  characteristics.  At  all 
frequencies  the  magnitude  difference  should  be  as  large  as  possible 
between  the  upper  and  lower  limits.  If  the  lower  limit  has  a  greater 
pole  to  zero  ratio  than  the  upper  limit,  the  magnitude  difference 
approaches  infinity  as  the  frequency,  <o,  approaches  infinity. 


Fig.  2-3  Pole-Zero  Pattern  for  a  Third  Order  Control  Ratio 
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A  Bode  Plot  (magnitude  plot  of  the  frequency  response)  for  each 
control  ratio  representing  a  time  response  bound  is  drawn.  Ihese  plots 
represent  the  system  specifications  in  the  frequency  domain  (Fig.  2-4). 
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Fig.  2-4  Design  Limits  for  a  Step  Input  (Frequency  Domain) 


These  equivalent  frequency  domain  specifications  generate  the  bounds  on 
the  loop  transmission,  L(jin). 
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The  system  specification  models  affect  the  final  design  in  several 
ways.  If  the  entire  set  of  acceptable  responses  is  not  considered,  the 
system  may  be  overdesigned,  thus  the  bandwidth  of  the  system  is  larger 
than  necessary.  If  the  frequencies  of  interest  change,  then  the  models 
need  changing,  else  the  design  may  be  overdesigned  or  not  acceptable. 


5.  Plant  Templates 

A  plant  template  is  a  plot  representing  the  range  of  uncertainty  in 
the  plant,  P,  at  a  given  frequency  (16:290).  For  example,  let  P(s)  - 
K/ 8 (s  +  a)  where  K  and  a  are  uncertain.  Given  2.0  <  K  <  8.0  db  and 
0.5  <  a  <  2.0  db  ,  an  infinite  number  of  possible  P's  exist  due  to  the 
uncertainties  in  K  and  a.  Generate  the  plant  template  by  plotting 
La[P(j<o)]  vs  AngpHjco)]  for  all  possible  values  of  R  and  a  at  a  suffi¬ 
cient  number  of  frequencies  in  the  range  of  interest.  This  generates  a 
-et  of  complex  numbers  with  finite  boundaries  for  any  frequency.  The 
templates  are  drawn  on  the  Nichols  Chart  because  the  chart  is  calibrated 
in  db  vs  please  angles  and  this  method  uses  the  Nichols  Chart.  Only  the 
boundary  of  the  template  need  calculating.  Find  the  boundary  of  the 
template  by  holding  all  but  one  variable  constant,  i.e.,  in  the  example 
above  set  K  ■  2  and  vary  a  incrementally  from  0.5  thru  2.0.  The  size 
of  the  increments  are  determined  by  the  designer.  If  the  shape  is 
simple,  use  fewer  points  (larger  increments)  or  if  the  shape  is 
irregular  use  more  points  (smaller  increments).  The  frequency  response 
at  iii  •  1  rad/sec  for  the  P's  obtained  above  provides  a  set  of  points 
from  A  (K  ■  2,  a  ■  0.5)  to  D  (K  ■  2 ,  a  ■  2)  on  the  Nichols  Chart 
(Fig.  2-5).  The  process  continues,  holding  all  but  one  variable 
coustant,  and  varying  that  one  variable  over  its  entire  range  of 
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Fig.  2-5  Plant  Templates  on  the  Nichols  Chart  (3:11-9) 
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uncertainty.  Increment  the  constants  while  the  changing  variable  is 
exercised  over  its  entire  range  of  uncertainty  for  each  set  of 
constants.  The  6ets  of  points  at  each  frequency  define  the  plant 
template  for  that  frequency.  Do  this  at  every  frequency  for  which  a 
plant  template  is  desired.  For  example,  given  a  -  0.5  ,  vary  K  from  2 
thru  8  to  obtain  Line  A  (a  *0.5,  K  ■  2)  to  B  (a  ■  0.5,  K*»8).  Now  set 
a  -  2.0  ,  and  vary  K  from  2  thru  8.  Complete  the  template  by  connecting 
the  four  corners  with  appropriately  curved  lines.  Make  a  plant 
template  for  many  frequencies,  the  designer  decides  on  the  intervals 
between  and  the  number  of  frequencies.  One  procedure  is  to  make 
templates  at  regular  intervals,  such  as  every  octave.  A  set  of 
templates  for  the  example  is  in  Fig.  2-5.  Note  the  size  changes  of  the 
plant  templates  (the  plant  template  for  0)  -  8.0  rad/sec  is  smaller 
than  the  plant  template  (o  *  1.0  rad/sec)  in  Fig.  2-5  at  the  different 
frequencies,  this  is  caused  by  the  range  of  uncertainty  in  P. 

6.  Nichols  Chart 

The  design  of  the  loop  transmission,  L,  and  the  prefilter,  F,  is 
performed  on  the  Nichols  Chart.  Let  F  be  equal  to  one  and  then  the  open 
loop  transmission  of  a  unity  feedback  system,  L  *  GP  ,  plots  onto  the 
horizontal  (magnitude  in  decibels)  and  vertical  (phase  angle  in  degrees) 
scales  of  the  Nichols  Chart.  At  any  given  frequency,  read  the  magnitude 
(in  decibels)  and  phase  angle  (in  degrees)  of  T  ■>  L/(l  +  L)  directly 
from  the  superimposed  curved  loci  [L/(l  +  L) ] .  Also  any  point  on  the 
loci  corresponding  to  the  magnitude  and  phase  angle  of  T  corresponds  to 
the  magnitude  and  phase  angle  of  L  on  the  horizontal  and  vertical  scales 
(5:332-334).  The  correspondence  between  L  and  T  is  important  because 

the  design  is  performed  using  this  relationship. 

16 


The  Nichols  Chart  is  also  used  for  disturbance  rejection.  Given 
that  «  1/(1  +  L)  (Eq.  2-2)  and  using  the  transformation,  L  «  1/m 
(1:152-155),  the  system  control  ratio  becomes  -  m/(l  +  m)  and  this 
is  of  the  same  form  as  ■  L/(l  +  L)  .  This  enables  the  design  of  the 
loop  transmission  inverse,  m,  directly  on  the  Nichols  Chart.  Turn  the 
Nichols  Chart  upside  down,  reflect  the  phase  angle  of  the  L  lines  about 
the  -180  degree  line  (i.e.,  -190  becomes  -170,  -200  becomes  -160,  etc.), 
and  reverse  the  signs  on  the  magnitude  lines  and  this  enables  the  design 
of  L  directly  on  the  Nichols  Chart.  The  horizontal  and  vertical  scales 
still  correspond  to  the  magnitude  and  phase  angle  of  L,  but  the 
superimposed  curved  loci  now  correspond  to  the  magnitude  and  phase  of 
1/(1  +  L)  (1:155).  For  the  disturbance  rejection  design,  only  the 
magnitude  is  necessary.  This  assumes  the  disturbance  has  its  maximum 
affect  possible  on  the  system  at  every  frequency.  Thus  ignore  the 
curved  lines  representing  the  phase  angles.  In  practice  for  a  disturbance 
rejection  design,  turn  the  Nichols  Chart  upside  down  and  modify  the 
scales  as  ^escribed  above.  Thus  design  L  directly  without  using  the 
dummy  variable,  m. 

7 .  Nominal  Plant 

Choose  an  arbitrary  specific  plant  parameter  set  to  generate 
VJ  *  a  reference  or  nominal  plant  used  in  the  design  method.  Pq 
generates  the  bounds  on  the  Nichols  Chart  for  the  nominal  loop  transmis¬ 
sion,  L  ,  where  L  ■  P  G  .  P  is  only  a  reference  and  therefore  is 
o  o  o  o 

anywhere,  even  outside  the  set  of  P's.  It  is  usually  convenient  to 
choose  a  P^  that  lies  within  the  set  of  possible  P's,  thus  being  a  point 
on  the  templates.  A  corner  of  the  template  is  convenient,  since  this 


transfer  function  is  already  known  and  calculated  for  each  template.  In 
the  above  example,  the  lower  left  hand  corner,  point  A,  is  chosen  as  the 


nominal  IPq  *  2/s(s  +  0.5)]  (Fig.  2-5).  It  helps  to  mark  the  reference 
point,  Pq,  on  every  template  because  the  Pq  chosen  may  not  remain  in  the 
same  relative  location  on  each  template. 

8.  Derivation  of  Bounds  for  the  Loop  Transmission,  L 

The  system  specifications  generate  the  bounds  for  each  loop  trans¬ 
mission,  L,  in  the  form  of  regions  permitted.  If  T(s)  uniquely  defines 
the  response,  y(t),  to  a  step  Input,  then  T(s)  is  completely  specified 
for  a  stable,  minimum  phase  (no  right  half  plane  poles  or  zeros)  system, 
by  the  magnitude  of  the  frequency  response,  |  T^( jco ) |  (14,16).  The 

frequency  response  of  |Y(j<u)|  varies  between  the  upper  (Ty)  and  the 
lower  (T^)  bounds.  For  the  example,  at  to  »  1  rad/sec  ,  assume  [_Y(jl)| 
varies  from  0.7  to  -0.8  db.  The  relative  variation  is  therefore 
I  (0.7)  -  (-0.8)|  *  1.5  db  .  The  allowable  relative  change  in  Y(jw)  at 
any  frequency  is  given  by 

ALm  [Kjw^]  -  ALm  [T^jc^)]  -  ALm  [T^jc^)]  (2-6) 

where  Ty(jto)  and  T^(jw)  are  the  frequency  domain  bounds  on  Y^(jco) . 

The  relationship  between  the  relative  change  and  the  control  ratio 
is  as  follows.  From  Fig.  2-1  and  Eq.  2-1,  ALm  (Y)  -  ALm  (.T)  = 

ALm  [FL/(,1  ♦  L)  ]  where  L  =  G  P  and  G  and  F  are  known,  i.e.  contain  no 
uncertainty.  Then; 


(2-7) 


ALm  [YCjoi^)]  «  ALm  [Tljo^)]  «=  ALm  C  L.C j / 1 1  +  LXjo^)}] 


Assuming  G  and  F  are  known  with  certainty  or  the  uncertainty  is 
negligible  (else  include  it  in  P) ,  the  relative  change  in  l-(ja)  i>  iS 
equal  to  the  relative  change  in  the  plant,  P^(jco^). 

ALm  [Ujw.)  ]  «=  ALm  [PCjoj^)]  (2-8) 

Parameter  uncertainty  in  P  is  the  cause  of  the  variation.  The  problem 
is  to  find  Lq  such  that  T(joj)  is  within  the  acceptable  range.  Thereby 
the  relative  change  requirements  on  the  closed-loop  response  are 
satisfied  over  the  entire  range  of  P's  uncertainty.  The  system 
specifications  are  the  requirements  on  the  closed-loop  response,  Y_ ( jto) , 
and  therefore  also  ]T(j<o)  (Eq.  2-7).  Next  transfer  the  constraints  for 

the  loop  transmission,  L(jo>) ,  to  the  Nichols  Chart  (14:18,16:291). 

9.  L(jai)  Bounds  on  the  Nichols  Chart 

The  relative  uncertainty  in  L  is  equivalent  to  the  range  of  uncer¬ 
tainty  in  P  above  (Eq.  2-8).  Remember,  the  plant  template  is  a  plot  of 
the  range  of  uncertainty  in  P  at  a  given  frequency.  Since  Lm  (L)  * 

Lm  (P)  +  Lm  (G)  and  Ang  (L)  *  Ang  (P)  +  Ang  (G)  ,  a  template  can  be 

translated  (but  not  rotated)  horizontally  or  vertically  on  the  Nichols 
Chart.  The  horizontal  and  vertical  translations  correspond  to  the  angle 
and  magnitude  requirements  on  G(jw)  at  a  given  frequency  (16:290). 
Drawing  lines  on  the  template  to  reference  the  horizontal  and  vertical 
helps  maintain  the  correct  orientation  on  the  Nichols  Chart  (see  Fig. 
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Using  the  templates  generated  for  the  example  in  the  previous 
sections,  place  the  template  for  oi  “  1  rad/sec  at  position  A  in  Fig. 

2-6.  The  area  covered  by  the  template  corresponds  to  the  variation  in  L 
due  to  the  uncertainty  in  P.  Recall  the  correspondence  of  L  and  T 
explained  in  Section  11-6.  Read  the  maximum  and  minimum  values  of  T 
covered  by  the  template  from  the  superimposed  curved  magnitude  contours 
of  the  Nichols  Chart  (the  contours  of  constant  Lm  [T(jaj)]).  If  this 
difference  (approximately  10  db)  between  the  maximum  and  minimum  values 
is  greater  than  the  allowable  difference  (1.5  db)  in  T  at  co  *  1  , 

ALm  [T(ja)^)]  (Eq.  2-7  and  Fig.  2-6),  shift  the  template  vertically 
upwards,  as  in  Fig.  2-6,  until  the  variation  covered  by  the  template  is 
equal  to  ALm  (T(j6),)]  (position  B  in  Fig.  2-6).  If  the  difference  is 
less  than  allowed,  then  translate  the  template  downward  until  the 
variation  covered  by  the  template  is  equal  to  the  difference  allowed. 
When  the  position  where  the  allowable  variation  and  that  covered  by  the 
template  are  equal,  transfer  the  point  for  the  nominal  plant  from  the 
template  onto  the  Nichols  Chart.  This  point,  read  from  the  horizontal 
and  vertical  scales  of  the  Nichols  Chart,  corresponds  to  a  bound,  B^(joj) , 
for  the  magnitude  and  phase  angle  values  of  the  L^.  The  nominal  loop 
transmission,  Lo(jw^),  is  given  by 

Lq (j^K)  -  G(jaK  )  Po(jco1)  (2-9) 

Repeat  this  process  horizontally  across  the  chart  for  different  values 

of  Ang(L  ).  The  points  form  a  curve,  B  ( jrj  ) ,  representing  the  boundary 
o  K  i 


lies  above  or  to  the  right  of  the  curve,  ),  the  variation  in  T  is 

less  than  or  equal  to  the  relative  change  allowed  by  the  system  specifi¬ 
cations  at  ok.  Repeat  this  process  at  every  frequency  of  interest, 
thereby  generating  a  set  of  design  bounds,  for  L  (jco) 

(16:291-292). 

In  order  to  effectively  reject  the  disturbance,  the  following 
inequality  must  be  satisfied: 


l/|l  +  L(jo)  |  <  [c(jeo)  | 


(2-10) 


where  |c(j<D)|  is  the  magnitude  of  the  boundary,  T^,  in  Fig.  2-4. 
Convert  the  magnitudes  to  decibels,  rearrange  the  terms  and  the 
inequality  becomes: 


Lm  [1  +  L(jco)]  >_-Lm  [C  ( jo>)  ]  (2-11) 

Place  a  template  on  the  modified  Nichols  Chart  (remember  the  Nichols 
Chart  is  modified  for  the  disturbance  rejection  case.  Section  11-6)  such 
that  its  lowest  point  rests  on  the  contour  of  constant  Lm  [1  +  L  (j  'A ) j 
equal  to  -Lm  [C(J<o^)]  ,  the  disturbance  bound  (B^(jW^)),  at  the 
frequency,  Transfer  the  nominal  point  to  the  Nichols  Chart  as 

before  and  slide  the  template  along  the  bound,  BD(jc^)  for  LQ(jw^). 
Repeat  this  process,  forming  a  set  of  bounds,  Bp(jco),  for  each  frequency 
of  Interest. 

Turn  the  Nichols  Chart  upright  again.  Currently  two  bounds  exist 

for  every  frequency  oi  interest.  The  command  response  bounds,  B  ,  and 
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the  disturbance  rejection  response  bounds,  B^  (see  Fig.  2-7).  Form  the 
composite  bound,  Bq,  by  eliminating  the  lower  of  the  two  bound  lines  at 
each  phase.  Thus,  the  composite  bound  consists  of  the  worse  of  the  two 
bounds  for  every  0)^.  Repeat  this  process  for  every  frequency,  thus 
forming  a  set  of  uominal  bounds,  BQ(jco)  (see  Fig.  2-8).  Use  these  worst 
bounds  to  shape  the  L^. 
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10.  Universal  high  Frequency  (UHF)  Bound 

The  Universal  High  Frequency  (UHF)  Bound  ensures  a  positive  gain 
and  phase  margin  for  the  loop  transmission,  L.  As  the  frequency 
approaches  infinicy,  the  plant  templates  become  narrower,  until  they  are 
considered  a  vertical  line.  The  allowable  variation  in  T  also  increases 
as  frequency  increases.  The  result  is  the  bound  of  tends  at  high 

frequency  to  a  very  narrow  region  around  0  db  and  -180  degree  point  (the 
origin)  on  the  Nichols  Chart.  The  UHF  Bound  is  used  to  avoid  placing  a 
closed-loop  pole  near  the  joj  axis,  which  results  in  a  oscillatory 
disturbance  response.  As  the  frequency  increases,  the  bounds  tend  to 
follow  the  ovals  encircling  the  origin.  Choose  one  of  the  ovals  near 
the  origin  as  the  bound.  For  this  example  choose  the  contour  of  constant 
magnitude  corresponding  to  5  db  (Fig.  2-8).  Find  the  template  with  the 
greatest  vertical  displacement,  Av  db ,  from  the  templates  oi  highest 
frequencies,  Av  is  also  determined  accurately  by  finding  the  maximum 
change  in  the  limit  of  Lm  [_P ( j  L>)  J  as  < .t  approaches  infinity.  Translate 
the  5  db  oval  down  the  length  of  the  template,  Av  db,  as  in  Fig.  2-8. 
This  defines  the  UHF  bound  (.14:20-22). 
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11. 

The  shaping  of  the  nominal  loop  transmission  to  satisfy  the 

boundaries  of  Lq  is  an  important  step  in  the  design.  Design  the  nominal 

L,  Lo>  to  satisfy  the  bounds  by  synthesizing  a  transfer  function  such 

that  the  magnitude  and  phase  does  not  violate  any  of  the  boundaries.  Do 

this  by  trial  and  error.  One  way  to  start  the  Lo  design,  is  to  use  P^. 

This  avoids  any  implicit  cancellation  of  roots  in  determining  G.  This 

also  simplifies  some  of  the  work  later,  since  L  ■  G  P  ,  where  P  and 

o  o  o 

L  are  known, 
o 

The  optimal  design  of  Lq  has  the  lowest  achievable  bandwidth  to 
solve  the  bounds.  The  optimal  Lq  has  the  value  of  the  bound  at  every 
frequency.  In  practical  designs,  the  goal  is  to  have  the  value  of  be 
higher  than,  but  as  close  as  possible  to,  its  corresponding  bound. 

Thus,  Lq  has  greater  bandwidth  than  the  optimal  Lq.  Therefore,  the 
design  is  a  tradeoff  between  getting  the  optimum  path  for  Lq  (minimum 
bandwidth)  which  requires  many  poles  and  zeros  and  using  a  simpler 
(larger  bandwidth)  with  fewer  poles  and  zeros. 

At  high  frequency,  relative  to  the  frequencies  of  interest,  add 
poles  (normally  complex  pairs)  to  lessen  the  bandwidth  of  the  design  as 
rapidly  as  possible. 

Figure  2-8  6hows  a  practical  design  for  the  above  example.  Any 
right  half  plane  (rhp)  poles  and/or  zeros  of  Pq  are  included  in  the  Lo 
to  avoid  any  attempt  to  cancel  the  rhp  poles  or  zeros  with  poles  or 
zeros  in  G.  For  another  example  and  a  discussion  of  loop  transmission 
shaping  see  Appendix  A. 


Shaping  of  the  Nominal  Loop  Transmission,  L 


Find  the  compensator  G  from  Lq. 


L  -  G  P 
o  o 


G  *  L  /P 
o  o 


(2-12) 


If  the  Lq  does  not  contain  the  roots  of  Po>  then  the  compensator  G  must 
cancel  them.  This  cancellation  occurs  only  for  the  purpose  of  design. 

In  actual  implementation,  cancellation  does  not  result  (nor  is  it 
necessary)  since  P  varies  over  the  entire  range  of  uncertainty.  With  G 
known,  the  loop  transmission,  L,  can  be  determined  for  any  particular 
plant  P(jO)). 

The  Lq  is  properly  shaped  when  it  meets  all  of  the  design  require¬ 
ments,  either  on  or  above  the  bound,  ®Q(jaii)»  at  each  frequency.  With  a 
properly  shaped  nominal  loop  transmission,  the  variation  in  T  (resulting 
from  the  uncertainty  in  P)  is  less  than  or  equal  to  the  allowable 
relative  change  in  T,  defined  by  the  system  specifications  (16:291). 

The  final  step  is  to  design  the  prefilter,  F. 

13.  Design  of  the  Prefilter,  F 

A  properly  designed  Lq  only  guarantees  that  the  variation  in 
) T ( j oo)  |  is  less  than  or  equal  to  the  allowed  variation.  The  prefilter 
positions  the  Lm  IT(jO))]  within  the  frequency  domain  system  specifications 
envelope.  For  the  example  the  magnitude  of  the  frequency  response  must 
remain  within  the  bounds  T  and  T  (see  Fig.  2-4).  The  bounds  are 

U  L 

redrawn  in  Fig.  2-9.  One  method  to  find  the  prefilter  is  as  follows. 

Place  the  nominal  point  of  the  oj  =  1  rad/sec  plant  template  on  the 
Nichols  Chart  where  the  L  (jl)  point  occurs.  Record  the  maximum  and 
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to  fit  within  the  bounds  and  T^.  For  example,  at  ^  =  1  rad/sec,  the 
actual  Lm  (T)  must  be  0.7  <  Lm  LT(jl))  <  -0.8  db.  From  the  plot  of  Lq, 
the  range  of  Lm  (T)  is  1.2  <  Lm  [T(jl)]  <  1.0  db.  To  lower  Lm  [T(jl>] 
into  the  desired  range,  the  prefilter,  Lm  (F)  required  is 


(0.7  -  1.2)  >  Lm  [F(j 1) ]  >  (-0.8  -  1.0) 


or  -0.5  >  Lm  IF(jl)]  >  -1.8  (see  Fig.  2-S)  (2-13) 


14.  Summary 

This  chapter  presented  an  overview  of  the  Quantitative  Feedback 
Theory  (QFT)  technique  of  Dr.  Horowitz  lor  a  S1S0  system.  The  technique 
is  performed  completely  in  the  frequency  domain,  requiring  only  simple 


Repeat  the  process  for  each  frequency  corresponding  to  the  plant 

templates  used.  Thus  the  difference  between  the  T„  and  T  curves  and 
r  U  max 

the  difference  between  the  T^  and  T  curves  indicate  the  requirements 

for  F  as  a  function  of  frequency  (Fig.  2-9). 

The  bounds  on  F,  (Lm  (T^)  -  Lm  (Tmav) ] >  Lm  (F)  >  [Lm  (T^)  - 

(T  ,  )]  ,  are  plotted  as  a  function  of  frequency  (Fig.  2-10).  Using  the 
min 

straight  line  approximation,  determine  the  transfer  function,  F,  whose 
magnitude  lies  within  these  bounds.  This  transfer  function  is  the 
prefilter  (16:301). 

The  single  loop  design  is  now  complete.  The  system  response  is 
guaranteed  to  remain  within  the  bounds  of  the  system  performance 
specifications,  provided  the  uncertainty  in  P  stays  within  the  range 
predicted  by  the  model  (16:288). 
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Fig.  2-10  Bounds  on  the  Prefilter,  F 

algebraic  calculations.  Nichols  Charts  and  Bode  Plots  are  used  exten 
sively  and  much  of  the  technique  is  performed  graphically. 

The  problem  is  defined,  then  the  system  specifications  are 
translated  into  the  frequency  domain.  These  specifications  generate 
limits  or  boundaries  of  the  frequency  response  for  the  system  control 
ratio  and  the  loop  transmission.  Two  compensator  elements,  G  and  F, 
control  the  system  response  to  inputs  and  disturbances  such  that  the 
system  response  meets  the  system  performance  specifications. 
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Ill.  QUANTITATIVE  FEEDBACK  THEORY: 
MULTIPLE  INPUT-MULTIPLE  OUTPUT  (MIMO)  SYSTEM 


1.  Introduction 


The  design  for  a  MIMO  system  involves  the  design  of  equivalent  S1S0 
system  loops.  The  design  process  for  these  individual  loops  are  the 
same  as  the  design  of  a  SISO  system  described  in  Chapter  II. 

In  general ,  an  n  x  n  MIMO  system  can  be  represented  in  matrix 
notation  as 


Pu 


(3-1) 


where  2.  m  vector  of  plant  outputs 
_u  ■  vector  of  plant  inputs 

_P  ■  plant  matrix  of  transfer  functions  relating  u  to  ^ 


The  P  matrix  is  formed  from  either  the  system  state  space  matrix  representa¬ 
tion  or  frdm  the  system  linear  differential  equations. 

Professor  Horowitz  has  shown,  by  using  fixed  point  theory  (15), 


that  the  inverse  of  the  P  matrix,  referred  to  as  contains  elements 


which  are  the  inverse  of  n  single  loop  transfer  functions  equivalent  to 


the  original  MIMO  plant.  Thus,  the  MIMO  problem  is  separated  into  n 


equivalent  SISO  systems  and  n  prefilter/disturbance  problems,  which  are 
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Consider  the  multiple  input-multiple  output  plant  of  Fig.  3-1.  The 
n  x  1  input  vector,  u_,  produces  an  n  x  1  output  vector,  the  n  x  n 
plant  matrix,  describes  the  relationship  between  and  u.  £  is  only  an 


Fig.  3-1  A  M1M0  Plant 

element  in  _P,  the  set  of  all  possible  plant  matrices.  Assume  the 
uncertainty  in  _P  can  be  determined,  or  at  least  bounded.  Also  assume 
the  input  and  output  vectors  are  of  the  same  dimension.  This  appears 
restrictive,  but  with  n  inputs,  at  most  only  n  outputs  are  independently 
controlled  (17:530-536).  If  the  model  defines  an  unequal  number  of 
Inputs  and  outputs,  the  first  task  is  to  modify  the  model  such  that  the 
dimensions  of  the  input  and  output  are  the  same. 

The  state  space  representation  for  a  MIMO  system  is: 


The  block  diagram  representing  these  equations  is  in  Fig.  3-2. 


Fig.  3-2  State  Space  Block  Diagram  for  Eq.  3-2 


If  the  plant  model  is  in  terms  of  n  coupled,  linear,  time  invariant 
differential  equations,  the  general  plant  model  for  a  second  order 
system  is: 

2 

(,a)y1  +  (bs  +  cs)y2  =  (f)uL  +  (g)u2 

(ds)yL  +  (e)y2  =  (h)u1  +  (i)u2  (3-3) 

where  a  thru  i  are  constant  coefficients,  the  y's  are  the  outputs  and 
the  u's  are  the  inputs.  In  matrix  notation  the  system  looks  like: 


bs  +  cs 


f  8 


h  i 


(3-4) 


Let  the  matrix  premultiplying  the  output  vector,  y_,  be  M  and  the  matrix 
premultiplying  the  input  vector,  11,  be  N..  The  system  now  looks  like: 


•/.  v.v. 


■>  V  V1 


V*3 

■  -  V. .•  «  V 


My  =  Nu 


(3-5) 
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The  plant  matrix  needed  is  defined  by: 


>  v  v 


Thus,  the  plant  matrix,  P^,  is 


Z  “  Pu 


(3-6) 


P  =  M_1N 


(3-7) 
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Another  method  used  to  generate  P^,  is  to  use  th.  state  space 
matrices  directly.  Although  any  number  of  states  are  represented, 
assume  that  the  input  and  output  vectors,  ii  and  are  of  the  same 

dimensions.  Assuming  the  system  is  linearized  and  the  A,  IJ  and  C 
matrices  are  time  invariant,  the  plant  matrix  is: 
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P  =  C[sl  -  A]  B 


(3-8) 


This  plant  matrix  is  actually  a  member  of  a  set  of  possible  plant 
matrices,  whose  set  is  caused  by  the  uncertainty  in  the  plant  parameters. 
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In  practice,  only  a  finite  set  of  £  matrices  are  formed  representing  the 
plant  under  varying  conditions. 


3.  MIMO  Compensation 

The  basic  MIMO  compensation  structure  for  a  two-by-two  MIMO  system 
is  in  Fig.  3-3.  The  structure  for  a  three-by-three  MIMO  system  is  in 
Fig.  3-5.  These  systems  are  similar  to  the  S1S0  system  of  Chapter  II. 

They  consist  of  the  uncertain  plant  matrix,  P ,  the  diagonal  compensation 
matrix,  (2,  and  the  prefilter  matrix,  F.  This  thesis  only  considers  a 
diagonal  £  matrix,  though  a  non-diagonal  G  matrix  gives  the  designer 
much  more  flexibility  in  his  design  (12:14-15).  A  more  detailed  schematic 
of  the  two-by-two  MIMO  system  is  given  in  Fig.  3-4  and  the  three-by-three 
MIMO  schematic  is  in  Fig.  3-6.  The  £  and  £  functions  are  the  same  as 
the  G  and  F  of  the  SISO  system. 

In  general  G  *  diag  £g±]  ,  JF  *  [f^]  and  P  *  Ip^]: 
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Fig.  3-5  MIMO  Control  Structure 
(Three-By-Three  System) 


Using  Eq.  3-9,  the  matrices  for  a  two-by-two  system  are 
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4.  Constraints  on  the  Plant  Matrix 

The  set  of  J?  must  be  tested  to  ensure  that  two  critical  conditions 
are  met  (15:86-90): 

1.  must  not  be  singular  for  any  combination  of  possible  plant 
parameters.  Thus  P  ^  exists. 

2.  As  s  +  lpllp22^  >  ^P12P2lJ  *or  Possible  plants  for  the 
two-by-two  case.  As  s  -*■  «,  IpuP22p33!  >  !pllp23P32^  + 

I p12p21p33 I  +  I P12P23P31 I  +  I p13p22p3ll  +  lp13p21p32l  for  <lLL 
possible  plants  for  the  three-by-three  case.  See  Reference  15 

for  higher  order  plants.  This  condition  is  necessary  only  if 

Eq.  3-7  and  the  original  method  (see  p.  42-43)  are  used  to 

generate  the  plant. 

The  first  condition  ensures  controllability  of  the  plant.  The 
inverse  of  P.  produces  the  effective  transfer  functions  used  in  the 
design.  If  the  £  matrix,  in  its  original  ordering  of  the  input  and 
output  matrices,  does  not  satisfy  the  second  condition,  then  a  reordering 
of  the  input  and  output  matrices  might  satisfy  the  condition. 

5.  Effective  SISO  Loops 

Define  the  matrix  Q*  =  P  '  whose  elements  are  q. 


2 1  n22  M2n 


P_1  *  o' 


(3-11) 
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Then  the  n  effective  transfer  functions  are  q  *=  1/q^,.  *  symbolically 
represented  by  5  *  1/5.  • 
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Reference  15  contains  the  derivations  and  proof  of  this  equivalence. 

2 

The  n  x  n  M1M0  system  is  transformed  into  n  SISO  problems.  The  general 

2 

transformation  result  of  n  SISO  system  loops  is  in  Fig.  3-7.  Figure 
3-8  shows  the  four  effective  SISO  loops  (in  the  boxed  area)  resulting 
from  a  two-by-two  system  and  the  nine  effective  SISO  loops  (13:682) 
resulting  from  a  three-by-three  system. 

Each  SISO  loop  in  Figs.  3-7  and  3-6  is  considered  as  an  individual 
SISO  design  problem,  which  is  solved  using  the  procedures  explained  in 


Fig.  3-7  Effective  SISO  Loops  (in  general) 


Chapter  II.  The  f 's  and  g’s  are  the  compensator  elements  of  the  J?  and  G 
matrices  described  previously. 

The  disturbances,  d^ ,  represent  the  interaction  between  the  loops: 
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iJ  i, 
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(3-15) 


where  the  numerator,  b^  ,  is  the  upper  response  bound,  (T^  or  in  Fig. 
2-4),  for  the  respective  input/output  relationship.  These  are  obtained 
from  the  design  specifications  ( 13 : 681-6b7) .  The  first  subscript,  k. 


refers  to  the  output  variable,  and  the  second  subscript,  j,  refers  to 


the  input  variable.  Therefore,  b.  .  is  a  function  of  the  response 

requirements  on  the  output,  y,  ,  due  to  the  input,  r..  The  lower  bound, 

J 

a^ ,  needs  defining  only  when  there  is  a  command  input. 
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For  the  disturbance  rejection  problem,  the  responses  must  be  less 
than  some  bound  (i.e.  £b^).  Thus  the  loop  equations  become: 
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where  -  L^.  Reorganizing: 
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using  Eq.  3-15,  and  rewriting  Eq.  3-17: 
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For  example,  in  a  three-by-three  system  for  the  first  loop,  L^, 
/ 

i  ■  1,  j  ■  2,  and  in  the  first  term  k  ■  2  and  in  the  second  term 
k  ■  3  ,  Eq.  3-18  becomes: 
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Remember  only  use  the  magnitude  in  the  disturbance  calculations.  This 
assumes  the  worst  case. 

Recently,  this  technique  was  improved.  This  improvement  involves 
modifying  the  q's,  using  the  information  from  the  previously  designed 
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loops,  the  g’s.  This  reduces  the  overdesign  inherent  in  the  early  part 
of  the  design  process.  The  final  loop,  is  designed  using  the  exact 
equation  (10:977).  The  exact  equation  represents  the  loop  and  the 
Interactions  caused  by  the  other  loops.  This  thesis  uses  the  original 
method  for  the  first  two  loops  and  the  improved  method  for  the  last 
loop. 

The  order  in  which  the  loops  are  designed  may  be  significant.  Any 
order  can  be  used,  but  some  orders  produce  less  overdesign  (less 
bandwidth)  than  others.  The  last  loop  designed  has  the  least  amount  of 
overdesign,  therefore  do  the  most  constrained  loop  first,  and  so  on 
through  all  of  the  loops. 

6.  Basically  Non-Interacting  (BNIC)  Loops 

A  basically  non-interacting  (BNIC)  loop  (13:679)  is  one  in  which 
the  output,  y^,  due  to  the  input,  r^>  is  ideally  zero.  Plant  uncertainty 
and  loop  Interaction  makes  the  ideal  response  unachievable.  Thus,  the 

system  performance  specifications  describe  a  range  of  acceptable  responses 

/ 

for  the  commanded  output  and  a  maximum  tolerable  response  for  the 
uncommanded  outputs.  The  uncommanded  outputs  are  treated  as  disturbances. 

7.  Summary 

This  chapter  describes  the  multiple  input-multiple  output  system 
and  the  plant  matrix.  Guidelines  for  finding  the  £  matrix,  which 
relates  the  input  vector  to  the  output  vector  are  given. 

The  separation  of  the  M1M0  system  into  several  S1S0  systems  is 
presented  using  the  £  ^ .  The  equivalent  S1S0  loops  are  designed 
according  to  the  S1S0  design  method  in  Chapter  II. 


IV .  AIRCRAFT  MODEL 


1.  Introduction 

The  aircraft  model  In  this  thesis  Is  a  hypothetical  aircraft 
designed  by  Lockheed,  Pratt  and  Whitney  and  Honeywell.  The  aircraft  was 
designed  as  part  of  the  study,  the  Flight  Propulsion  Control  Coupling 
(FPCC)  and  Dynamic  Interaction  Investigation  (7,8).  This  hypothetical, 
or  "paper"  aircraft  is  referred  to  as  the  FPCC  aircraft  (Fig.  4-1). 

The  aircraft's  basic  design  criterion  is  for  a  future  multimission 
fighter  whose  primary  function  is  supersonic/transonic  air  superiority. 
Its  secondary  function  is  transonic  close  air  support.  The  aircraft's 
dash,  or  top  speed  is  Mach  2.2  to  3.  The  aircraft's  takeoff  weight  is 
34,000  pounds  and  the  maneuvering  weight  is  29,000  pounds.  The 
maneuvering  wing  loading  is  44.3  pounds  per  square  foot.  It  is  capable 
of  conventional  or  short  takeoffs  and  landings  and  has  relaxed  static 
stability  and  direct  side  force  control. 

The  aircraft  uses  many  flight  control  surfaces  to  meet  the  require¬ 
ments.  The  primary  flight  control  surfaces  are  maneuvering  flaps, 
leading  edge  flaps.  Jet  flaps,  canards,  canard  flaps,  chin  canards, 
upper  rudder,  lower  rudder,  and  ailerons  IFig.  4-2).  The  vertical 
(chin)  canards,  rudder  and  ailerons  generate  the  forces  for  lateral 
control. 


2 .  Aircraft  Description 

This  thesis  involves  designing  a  lateral  flight  controller.  The 


equations  of  motion  are  not  explicitly  used  by  the  author  because  the 


computer  simulation  program,  FPCCSIM,  is  a\„ilable.  Appendix  A  oi 
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and  Function  (7:35) 


Reference  8,  pages  157-163,  contains  the  equations  of  motion.  FPCCS1M 
generates  the  state  space  model  of  the  aircraft  at  the  flight  conditions 
of  interest.  Note  the  aircraft  is  not  designed  for  decoupled  maneuvers, 
though  it  does  have  the  capability.  Also  the  aircraft  is  statically 
unstable  at  mach  0.9. 

The  FPCSSIM  computer  program  was  written  as  part  of  the  FPCC  study. 
The  FPCCSIM  simulation  uses  a  non-linear  6  degree-of-freedom 
perturbation  model.  This  is  the  same  program  used  by  previous  AFIT 
thesis  students  in  their  work  using  other  design  techniques. 

The  command  inputs  available  are  roll  angle  (4>)»  side  velocity  (v) 
and  yaw  rate  (r).  The  commanded  inputs  are  limited  in  this  thesis  to 
deflections  of  the  ailerons,  the  rudder  and  the  vertical  canard.  The 
outputs  measured  are  roll  angle,  roll  rate  (p),  side  velocity  and  yaw 
rate. 

As  stated  above,  the  aircraft  model  is  in  the  state  space  form: 

/  A  *  Ax  +  Bu  (4-1) 

X  =  Cx  (4-2) 

For  example,  the  Eq.  4-3  and  4-4  are  outputs  from  FPCCSIM  for  the  slow 
speed  condition  (mach  0.6  at  sea  level). 

The  state  vector  and  the  control  vector  are: 
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where  the  plant  state  space  elements  (x)  are: 

P 

body  axis  roll  rate  (rad/sec) 

r 

body  axis  yaw  rate  (rad/sec) 

V 

body  axis  side  velocity  (ft/sec) 

roll  angle  (rad) 

* 

yaw  angle  (rad) 

y 

cross  range  position  reference  to  initial  body  axis  (ft) 

q 

body  axis  pitch  rate  (rad/sec) 

g 

body  axis  pitch  rate  (rad/sec) 

w 

body  axis  vertical  velocity  (ft/sec) 

u 

body  axis  forward  velocity  (ft/sec) 

e 

pitch  angle  (rad) 

h 

altitude  (ft) 

X 

uown  range  position  re i  deuce  to  initiuj  hoc.  liU 
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and  the  control  input  elements  (jr)  are: 

6  sum  of  both  aileron  deflections  (deg) 

<5  horizontal  canard  deflection  (deg) 
c 

<5  maneuver  flap  deflection  (deg) 

6^  rudder  deflection  (deg) 

6  vertical  canard  deflection  (deg) 

vc 

6^  jet  flap  deflection  with  respect  to  fuselage  line  (deg) 
net  thrust  engine  1  (lbs) 
net  thrust  engine  2  (lbs) 

CDI^  inlet  drag  coefficient  1 

CDl^  inlet  drag  coefficient  2 

Appendix  B  contains  the  FPCCS1M  generated  matrices.  FPCCSIM  also 
gives  the  trim  condition  of  the  aircraft  for  the  flight  condition.  For 
this  flight  condition,  0.6  mach  at  sea  level,  the  angle  of  attack  is 
1.03  degrees  at  0.6  mach  with  a  dynamic  pressure  of  533.04  lbs/sq  in  and 
a  forward  velocity  of  670.5  ft/sec. 

3.  Lateral  Equation  Model 

The  above  matrices  are  simplified  by  noting  they  contain  both  the 

lateral  and  longitudinal  terms.  The  lateral  terms  in  the  plant  state 

space  matrix.  A,  are  p,  r,  v,  $,  y,  and  ip.  Both  y  and  tp  are  eliminated 

with  little  degradation  of  the  model.  The  lateral  terms  in  the  control 

input  matrix,  B,  are  6,6,6  ,  F.  ,  F,,,  CD1.  and  CDI„.  The  first 

—  a  r  vc  1  2  1  2 

three  terms  are  the  significant  contributors  in  the  lateral  mode.  Thus, 
the  model  is  simplified  to  a  4  x  4  plant  state  space  matrix  and  a  4  x  3 
control  input  matrix.  Simplifying  and  rewriting: 
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The  other  lateral  models  are  in  Appendix  B. 


4.  Flight  Conditions 

The  three  flight  conditions  used  represent  the  aircraft's 
operational  envelope.  There  is  a  slow  speed  condition.  Flight  Condition 
#1  (FC  #1),  which  is  at  maximum  weight  (34,000  pounds)  at  sea  level  with 
full  dry  thrust  (no  afterburner)  and  a  speed  of  0.6  mach.  The  cruise 
condition.  Flight  Condition  it 2  (FC  it2) ,  is  at  maneuvering  weight  (29,000 
pounds)  at  30,000  feet  and  a  speed  of  0.9  mach.  The  third  case,  Flight 
Condition  #3  (FC  it 3),  is  at  maneuvering  weight  at  40,000  feet  and  a 
speed  of  2.3  mach. 


The  maneuver  performed  is  a  horizontal  translation.  This  maneuver 
requires  displacing  the  aircraft  laterally  without  a  yaw  rate  or  bank 
angle.  The  maneuver  is  performed  by  commanding  a  side  velocity.  The 
roll  angle  and  yaw  rate,  generated  by  coupling  with  the  side  velocity, 
are  commanded  to  zero.  Any  roll  angle  or  yaw  rate  resulting  from  the 
coupling  effects  is  treated  as  a  disturbance,  since  it  is  desired  that 
both  remain  zero. 

The  three  flight  conditions  and  the  maneuver  are  the  6ame  that  a 
previous  thesis  (2)  used  with  a  different  design  technique.  Part  of 
this  thesis  compares  the  results  from  that  thesis  to  this  one  (Chapter 
Vlll).  Therefore,  the  same  state  space  plant  matrices  must  be  used. 
Thus,  the  flight  conditions  are  matched  such  that  the  three  plant  state 
space  matrices  are  the  same. 
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V.  DESIGN  OF  THE  FIRST  TWO  LOOPS 


1.  Introduction 

This  section  presents  in  detail  the  first  two  loops  designed  for 
the  flight  controller.  First  the  plant  matrices  are  generated  and  then 
the  system  specifications  are  chosen.  Next  the  order  in  which  to  design 
the  loops  is  selected  and  the  plant  templates  generated.  The  design 
phase  begins  by  generating  bounds  for  the  first  two  loops  and  then 
designing  the  open  loop  transfer  functions. 

2.  The  Problem 

Design  a  lateral  flight  controller  for  the  horizontal  translation 
maneuver  using  the  QFT  technique.  The  plant,  the  Flight  Propulsion 
Control  Coupled  (FPCC)  aircraft,  is  modeled  as  a  three  input-three 
output  system  in  the  lateral  axis  (Section  IV-3) .  From  Section  IV-4, 
the  only  input  is  side  velocity,  v.  Optimally,  the  other  two  control 
variables  whose  inputs  are  zero  (roll  angle,  4>  and  yaw  rate,  r)  remain 
zero.  The  system  specifications  are  derived  from  the  responses  of  a 
previous  thesis  (2). 

3.  The  Plant  Matrices 

The  QFT  technique  requires  the  plant  description  in  the  transfer 
function  form.  The  FPCCSIM  computer  program  generates  state  space 
models  for  each  of  the  flight  conditions  (FC).  Another  computer  program 
generates  the  transfer  function  plant  matrix  from  the  state  space 
matrices  for  each  FC. 
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Appendix  B  contains  the  full  aircraft  state  space  matrices.  Since 
this  design  includes  only  the  lateral  axis,  the  full  state  space 
matrices  are  unnecessary.  Simplify  the  full  state  space  matrices  as 
explained  in  Section  1V-3,  to  a  three-by-three  system.  These  sets  of 
matrices  are  the  plant  models.  Appendix  B  lists  the  simplified  lateral 
mode  state  space  matrices.  A  computer  program  generates  the  plant 
transfer  function  matrices,  £  (in  Appendix  B),  using  the  lateral  mode 
state  space  matrices  as  input.  With  the  plant  in  the  form  of  transfer 
function  matrices  (one  for  each  flight  condition  considered),  calculate 
the  (£' s  using: 


(3-13) 


a-  lqtj]  -  i/cf  L 1  -  U/qyl 


(3-14) 


Appendix  B  contains  the  (^'s  for  all  ol  the  flight  conditions. 


4.  System  Specifications 

The  system  specifications  on  t^  where  T  *  [t  ]  »  are  formulated 
in  the  frequency  domain.  This  thesis  uses  the  system  specifications 
from  the  response  functions  of  a  previous  thesis  (2).  These  responses 
are  used  because  a  military  specification  for  the  horizontal  translation 
maneuver  did  not  exist.  Based  on  the  response  curves  of  the  previous 
thesis  (Figs.  5-1  and  5-2),  synthesize  the  frequency  domain  system 
specifications  from  the  time  domain  responses  (see  Section  11-4) .  Table 
5-1  contains  the  specifications  in  the  time  domain.  The  b's  are  defined 
the  same  as  in  Section  111-5  (i.e.  b ^  defined  as  the  maximum  acceptable 
response  of  the  third  loop  due  to  an  input  at  the  second  loop) . 


Table  5-1  Time  Domain  System  Specifications 
(approximated  from  previous  thesis  (2)) 


t 

t 

t 

M 

F 

r 

P 

s 

P 

V 

22 

1.3 

3.7 

6.0 

1.33 

1.0 

32 

1.0 

1.0 

5.6 

0.03 

0.0 

12 

1.2 

1.2 

9.8 

0.05 

0.0 

22 

2.2 

3.8 

1.0 

1.0 

Table  5-2  Equivalent  Transfer  Function  Specifications 
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Using  these  time  domain  system  response  curves  of  the  previous 
thesis,  generate  the  equivalent  frequency  domain  system  specifications 
(using  the  method  in  Section  11-4).  Some  leeway  in  defining  the  system 
specifications  exists.  The  synthesis  of  the  equivalent  transfer  functions 
includes  compromises.  A  transfer  function  which  satisfies  all  of  the 
time  domain  specifications  requires  many  iterations  to  synthesize.  Thus 
choose  a  few  of  the  specification  criteria  as  important  ones  and  compare 
these  important  criteria  of  the  equivalent  transfer  functions  to  tile 
desired  time  domain  specifications.  The  important  specifications  chosen 
are  the  peak  value  (M^),  the  final  value  (F^)  and  the  rise  time  (t^). 

Using  these  three  specification  criteria,  the  equivalent  transfer 
functions  are  synthesized  faster  and  are  of  a  simpler  order  than  if  all 
of  the  time  domain  specifications  are  matched.  This  is  a  compromise 
between  the  work  and  time  necessary  to  generate  the  transfer  functions 
and  how  close  the  transfer  function  specifications  are  to  the  time 
domain  specifications. 

An  iterative  process  is  used  to  generate  the  equivalent  transfer 
functions.  As  an  example,  the  synthesis  of  the  loop  2  lower  bound,  a^* 
follows.  Since  only  specific  responses  are  given,  a  range  of  acceptable 
responses  is  generated  based  on  the  responses  and  experience.  Using  the 
criteria  from  the  previous  thesis  (Fig.  5-1  and  5-2),  a ^ ,  the  lower 
bound,  is  overdamped  with  no  overshoot.  The  response  starts  at  zero  and 
settles  to  a  final  value  of  one.  The  rise  time  is  approximately  2  sec 
and  the  settling  time  is  about  4  sec  (Table  5-1). 

Using  the  known  information  as  a  starting  point,  synthesize  a^* 

Since  it  is  overdamped,  try  a  damping  ratio  greater  than  one.  Since  the 
final  value  is  one,  the  steady  state  response  must  be  one.  The 


frequency  response  characteristics  of  the  final  transfer  function  is 

also  important.  Remember  the  response  should  decrease  as  fast  as 

possible  in  the  frequency  domain  (it  is  a  lower  bound)  thereby  giving 

the  largest  possible  region  between  the  upper  bound,  and  a ^ • 

Thus,  desire  a  transfer  function  with  three  more  poles  than  zeros,  or  at 

least  two  more  poles  than  zeros. 

Start  the  iterations  by  assuming  a  set  of  poles  with  a  natural 

frequency,  such  as  five,  and  a  damping  ratio  greater  than  one.  Keep  the 

natural  frequency  constant  and  vary  the  damping  ratio.  Start  with  a 

2 

damping  ratio  of  one.  Call  this  a  .,  a  ,  =  25/ (s  +  10s  +  25).  The 

xl  xl 

time  response  criteria  of  a  ,  are  H  *  1,  F  **  1,  t  «■  0.672  and  t  ■ 
r  xl  p  v  r  s 

1.167.  The  responses,  tf  and  tg  are  too  fast.  Try  a  damping  ratio  of 

two,  a  0  -  25/ (s2  +  20s  +  25).  This  yields  M  «  1,  F  -  1,  t  «  1.646 
xx  p  v  r 

and  t  *  2.9756,  where  t  and  t  are  still  too  fast,  but  they  are  closer 
8  r  s 

to  the  desired  response  than  a  Thus  increase  the  damping  ratio 

2 

again.  Try  a  damping  ratio  of  three,  a^  ■  25/ (s  +  30s  +  25).  This 

yields  M  «=  1,  F  *=  1 ,  t  *  2.562  and  t  ■  4.595,  where  t  and  t  are 
p  v  r  s  r  s 

now  too  slow.  Of  the  three  functions,  a  „  is  the  best. 

xl 

The  designer  decides  this  is  not  close  enough,  therefore  the 

Iteration  process  continues.  Try  adding  another  pole  to  a^*  ax2  * 

25/ (s2  +  20s  +  25) (s  +  1)  .  This  yields  M  -  1,  Fy  -  1,  tf  -  2.9499  and 

t  *  5.1998.  Now  t  and  t  are  too  slow.  Try  a  pole  at  two,  a  „  * 
s  r  s  j  v  x2 

50/ (s2  +  20s  +  25)  (s  +  2).  This  yields  M  *=  1,  F^  =  1,  t^  =  2.116  and 

t  “  3.76,  which  are  close  to  the  desired  specifications.  Try  a  pole  at 
8 

three,  a  „  *  75/(s2  +  20s  +  25) (s  +  3).  This  yields  M  =  1,  F  =1, 
x  z  p  V 

t  *=  1.8866  and  t  =  3.416,  where  t  and  t  are  too  last.  ihe  a  „ 
r  s  r  s  x2 
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function  with  the  pole  at  -2  gives  the  best  results.  The  designer 

O 

considers  the  time  response  from  this  function,  =  50/ (s^  +  20s  + 

25) (s  ♦  2)  close  enough.  Thus  the  iterative  process  ends  and  a 22  is 
synthesized. 

The  transfer  function  chosen  is  not  the  only  possible  solution. 

The  designer  decides  when  the  equivalent  transfer  function  specifications 
match  the  desired  specifications  well  enough.  The  harder  the  bounds  are 
to  satisfy,  the  more  careful  the  synthesis  of  these  equivalent  transfer 
functions  must  be.  The  QFT  technique  uses  the  values  of  the  equivalent 
transfer  functions  in  the  frequency  domain,  and  these  values  are  sometimes 
referred  to  as  the  equivalent  frequency  domain  system  specifications. 

Use  the  same  procedure  to  synthesize  equivalent  transfer  functions 
for  the  three  other  bounds.  The  upper  bound  for  L2,  b^,  also  settles 
to  a  final  value  of  one,  thus  it  too  has  a  steady  state  value  of  one. 

The  other  two  bounds,  b^  and  b^2>  settle  at  zero.  They  need  a  zero  (a 
numerator  root)  at  the  origin. 

Some  of  the  lesser  important  specification  criteria  are  violated. 
These  violations  are  minor,  less  than  5%  (Table  5-2)  and  believed  to 
only  slightly  degrade  the  model.  This  approximation  does  not  affect  the 
comparison  between  the  two  theses.  Three  decibels  (3  db  gain  margin)  is 
selected  as  the  Universal-High  Frequency  (UHF)  Bound. 

The  equivalent  transfer  functions  obtained  are: 


-TV  ■  j  .ji  .j,  .j,  U  ■  L  ■  l  ■  >  ■  L  W  W 11 W  A  l*  A  WH  'A  TTO  "•"■  '*•  ■  ^ 


-0.33  s  (s  -  l)(s  +  1.4) 
b32  =  (s  +  0.3) (s  +  1.5)  (s  +  2)  (s  +  4) 

.  0.36  s(s  -  2)  (s  +  0.25) 

b12  "  (s  +  0. 15)  (s  +  0.5)  (s  +  2)  (s  +  5) 

_ 50 _ 

a22  ”  (s  +  18.66) (s  +  2)(s  +  1.34) 

This  iterative  process  is  also  explained  in  Section  11-4.  Fig.  5-3 
shows  the  Bode  Plot  of  the  equivalent  frequency  domain  system  specifica¬ 
tions.  Note  the  upper  level  and  lower  bound,  b ^  and  a^*  define  an 
envelope  which  the  system  response  is  between.  Notice  both  of  the  other 
bounds  (which  are  effectively  upper  bounds)  are  below  the  upper  bound  of 
the  desired  system  response,  Also  the  two  bounds,  b^2  and  b^>  are 

below  0  db,  which  means  these  responses  are  attenuated.  The  larger  the 
difference  between  the  envelope  of  acceptable  responses  (the  region 
between  b^  and  a^)  and  the  bounds  (b^  aD(*  b12^  tbe  “desired 
responses,  the  easier  to  perform  the  design.  Since  the  QFT  technique 
uses  the  system  specifications  explicitly,  the  designer  knows  during  the 
design  phase,  when  a  specification  is  violated  and  can  take  appropriate 
action.  Thus  using  these  system  specifications  this  design  should  have 
at  least  similar  system  responses  as  the  previous  thesis  (2). 

5.  Choosin,,  Which  Loops  to  Design  First 

kith  the  system  specifications  and  the  plant  matrices  in  the 
required  form,  the  design  phase  begins.  Design  the  loop  with  the 
highest  bandwidth  requirements  last,  since  the  last  loop  designed 


contains  the  least  amount  of  oveiuesign  compared  to  the  uthcr  loops 


(Section  III-5).  Remember  the  maneuver  performed,  a  horizontal  trans¬ 
lation  has  only  one  non-zero  commanded  input,  the  side  velocity,  v. 

Loop  2  (L2)  includes  the  input  command  and  the  coupling  effects  from  the 
other  two  loops,  and  the  other  two  loops,  and  L^,  contain  only  the 
interaction  (or  coupling)  effects  as  input  (Section  III-5)  .  Thus,  of 
the  nine  possible  loops  in  a  three-by-three  system  (see  Fig.  3-8),  this 
system  requires  only  three  loops  designed  (Fig.  5-4).  Design  L2  last, 
since  this  loop  is  expected  to  require  more  bandwidth  than  the  other  two 
loops.  Also,  L2  requires  more  calculations  because  it  involves  both 
command  input  and  "disturbance"  input  terms,  whereas  the  other  two  loops 
only  contain  "disturbance"  terms.  Design  either  of  the  other  two  loops, 
(roll  angle,  <J>)  and  (yaw  rate,  r),  first  since  this  thesis  uses 
the  non-interacting  (the  original,  not  the  improved;  see  Section  III-5 
and  III-6)  method  for  the  first  two  loops.  Therefore  the  design  of 
these  two  loops  do  not  explicitly  affect  one  another. 

Since  both  loops,  L ^  and  L^  use  the  same  design  method,  this 
chapter  presents  only  one  loop,  L^.  Appendix  C  contains  the  loop 
design. 

6.  Plant  Templates 

Generate  the  plant  templates  for  each  frequency  of  interest. 

Choose  the  frequencies  of  interest  as  follows.  Flight  control  systems 
usually  operate  between  one  and  four  hertz,  which  is  6.28  to  25.08  rad/sec 
Therefore  use  a  decade  on  either  side.  This  choice  is  at  the  designer's 
discretion  based  upon  the  amount  of  uncertainty  at  the  frequencies  of 
interest  (Section  III-5).  Use  as  the  intervals,  octaves  between  the 
decades  (also  decided  by  the  designer) .  The  frequencies  of  interest 
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Fig.  5-4  Equivalent  SISO  Systems  for  a  3  x  3  M1M0  System 


chosen  are;  0.1,  0.2,  0.4,  0.8,  1.0,  2.0,  4.0,  8.0,  10.0,  20.0,  40.0, 
80.0  and  100.0  rad/sec.  Table  C-l  lists  the  values  of  (Table  C-2 
contains  q^ 3).  uses  the  q^  term  as  the  plant  (L^  uses  q^) . 

Fig.  C-l  shows  a  Bode  Plot  of  the  q^  and  q^  terms.  These  magnitudes 
generate  the  plant  templates  as  explained  in  Section  11-5.  Fig.  C-2 
shows  the  plant  templates  ([q^]  sets)  for  and  Fig.  C-3  for  L^. 

With  the  plant  templates  generated  and  the  system  bounds  defined. 


calculate  the  bounds  for  each  loop. 


Bounds  for  the  First  Two  Loods 


and  are  effectively  disturbance  rejection  loops  (no  command 
input).  The  system  equations  for  and  (Fig.  5-4  and  Eq.  3-20)  are: 


li  +  4 1  1  < 


qn 

.  .  qll 

b22 

qi2 

+  b-„ 

32 

(5-1) 


b  211  ^33 

,  b22  q  +  b12  q_ 

U  +  L3I  1  < - 1--32|  - LJll 


(5-2) 


These  expressions  require  four  ratios,  I  *  Iqll/fqi3l  * 

1*133/ *^32  I  an<*  lq33^q31  I  (tables  C-l  thru  C-6) .  Remember  from  Section 
II-8,  Iq^/q^l  is  a  ratio  of  magnitudes,  not  an  equation.  Only  use  the 
magnitude  of  the  q  terms.  For  example,  for  FC  if 3,  at  a  frequency  of  10 
rad/sec;  qu  -  -51.3  db,  q12  -  7.6  db  and  q  -  -38.3  db.  Thus: 


=  -58.9  db 


Oj  3  =  -13.0  db 


This  in  effect  considers  the  worst  case;  the  q  terms  have  their  maximum 
affect  on  the  nominal  bounds.  The  system  expressions  also  require  the 


magnitudes  of  the  bounds  b22,  ^32  an<*  ^12"  C-7  lists  the  magnitudes 

of  these  terms.  Substituting  these  values  into  Eqs.  5-1  and  5-2  yields 
the  bounds  (in  decibels)  for  the  appropriate  loop.  For  example,  for 
FC  #3,  at  oj  -  10  rad/sec,  b22  =  -17.4  db,  b32  «=  -30.4  db,  b  -  -29.9  db, 

| q  1 1  /q  1 2!  “  -58.9  db  andlq^/q^l  *  -13.0  db.  Substituting  these  values 
into  Eq.  5-1: 

1,  .  T  i-l  .  -17.4  (-58.9)  +  -30.4  (-13.0) 

|l+Lll  1  12979 


and  simplifying: 


,  ,  -76.3  +  -43.4  _ 

1  Ll'  - - =2979 


-43.2 

-29.9 


-13.3  db 


Table '5-3  lists  the  boundary  values  for  the  loop  (Table  C-8 
lists  the  loop).  Generate  the  bounds  for  the  first  loops  designed. 

As  explained  in  Section  II-8  use  the  highest  bound  of  the  three  flight 
conditions  at  each  frequency  as  the  nominal  bound.  Now  transfer  the 
nominal  bounds  for  to  the  Nichols  Chart,  as  explained  in  Section  11-8 
(Fig.  5-5). 

Remember  the  design  performed  is  actually  for  a  nominal  (reference) 
loop  transmission,  L^,  where  *=  P^q  (Section  11-11).  This  design 
method  uses  a  nominal  plant,  Pq  (Section  II-7)  which  can  be  anywhere. 

For  this  example,  P  is  chosen  as  flight  condition  (FC)  g 3  (Fig.  C-2). 


I  iLQJLiU 


r'l 


Table  5-3 


Loop  Bounds 


Nominal 


Freq. 

FC  #1 

FC  #2 

FC  # 3 

Bound 

*> 

(rad/sec) 

(db) 

(db) 

(db) 

(db) 

% 

0.1 

20.52 

23.52 

32.62 

32.62 

0.2 

16.21 

23.04 

31.70 

31.70 

& 

0.4 

10.66 

20.41 

28.76 

28.76 

0.8 

4.82 

15.30 

23.22 

28.22 

1.0 

3.08 

13.36 

21.07 

21.07 

:£ 

2.0 

-2.99 

5.78 

12.74 

12.74 

V*. 

4.0 

-10.76 

-4.32 

2.20 

2.20 

V 

"  hi 

8.0 

-18.88 

-15.56 

-9.57 

-9.57 

•  * 
> 

10.0 

-21.27 

-19.09 

-13.39 

-13.39 

20.0 

-28.02 

-28.93 

-24.72 

-24.72 

.V 

40.0 

-34.25 

-36.84 

-34.64 

-34.25 

80.0 

-35.23 

-43.50 

-42.66 

-35.23 

> 

100.0 

-42.28 

-45.53 

-44.92 

-42.28 

'■i; 

Thus ,  in 

this  example 

is 

the  same  as  i^ 

for  FC  #3  [(L1 

)31. 

This 

*  , 

S 

V 

tm* 

reference 

loop  transmission 

is  a  means  for  : 

finding  the  loop  compensator. 

•\* 

Gx  (Section  11-12).  With 

G1 

known  from 

•  L10  G1  P10' 

the  loop 

.y 

transmissions  for  any  plant 

can  be  derived. 

‘Vi  -  G1  (P1 

V 

Thus  in 

-• 

synthesizing  a  nominal  loop 

transmission  it 

is  the  compensator 

which  is 

being  designed. 

f- 

Use 

the  nominal  bounds 

for  designing  the  L^q.  The  nominal  bounds 

.v' 

V 

are  the  worst  bounds  from 

the  set  of  bounds 

generated  for 

the 

L^  loop. 

* 

(< 

From  Table  5-3,  the  bounds  for  L^,  at  each  of  the  three  flight  conditions 
of  interest  are  listed.  If  the  worst  bounds  are  not  used,  the  designed 
G  would  not  guarantee  the  loop  transmissions  to  satisfy  their  system 
specifications  at  each  flight  condition,  because  one  flight  condition 
might  violate  the  system  specifications.  Thus  use  the  nominal  bounds  to 
synthesize  a  and  thereby  design  such  that  the  system  specifica¬ 


tions  are  met  at  all  of  the  flight  conditions  of  interest. 


For  example,  transfer  the  <d  =  2  rad/sec  bound  onto  the  Nichols 
Chart.  Take  the  template  corresponding  to  tu  =  2  rad /sec  for  (Fig. 
C-2) .  Since  this  is  a  disturbance  rejection  problem,  invert  the  Nichols 
Chart  (Section  11-9)  and  the  template.  The  reference  point,  or  nominal 
plant  used  (Section  11-7),  is  FC  #3.  With  the  template  inverted,  the 
reference  point  is  at  the  bottom  of  the  template.  The  nominal  bound 
calculated  for  L^(j2.0)  is  12.7  db  (Table  5-3).  Put  the  lowest  point  of 
the  template  (in  this  example,  the  reference  point)  on  the  Nichols  Chart 
at  the  bound  (12.7  db  curve).  Remember  the  horizontal  and  vertical 
reference  on  the  template  must  stay  parallel  to  the  horizontal  and 
vertical  lines  of  the  Nichols  Chart.  Slide  the  template  along  the 
bound,  12.7  db,  transferring  the  reference  mark  from  the  template  to  the 
Nichols  Chart  at  enough  intervals  to  enable  drawing  a  representative 
curve  through  these  marks.  Since  the  lowest  point  of  the  template  is 
the  reference  point,  the  curve  drawn  is  the  12.7  db  curved  line  (Fig. 
5-5),  this  is  the  design  bound.  When  the  UHF  Bound  (the  3  decibel  curve 
chosen  in  Section  V-A) ,  becomes  the  limiting  bound  (for  L^,  when  0)  >  8 
rad/sec)  notice  the  template  causes  an  offset  from  the  UHF  Bound  (Fig. 
5-5).  If  any  other  nominal  plant  was  chosen  this  offsetting  effect 
would  occur  at  the  other  frequencies  as  well.  Note  from  Fig.  5-5  at  low 
frequency  the  bounds  tend  towards  a  straight  line  across  the  Nichols 
Chart,  but  tend  to  encircle  the  UHF  Bound  as  the  frequency  increases. 

The  UHF  Bound  keeps  the  open  loop  transfer  function  design  from  going  to 
the  origin,  -180  degrees  and  0  db,  which  causes  an  instability.  The 
design  bounds  are  in  Fig.  C-5. 


Do  this  at  all  frequencies  of  interest.  Every  frequency  cf  interest 


has  a  bound.  In  the  example,  as  the  frequency  increases,  the  nominal 
bounds  are  encompassed  by  the  UHF  Bound.  Therefore  use  the  UHF  Bound 
for  the  design  bound  at  these  frequencies  because  the  UHF  Bound  is  more 
restrictive. 

8.  Shaping  the  Open  Loop  Transfer  Function 

With  the  bounds  drawn,  shape  the  nominal  open  loop  transfer  function, 
L^q,  as  explained  in  Section  11-11  and  Appendix  A.  Iterations  used  to 
shape  are  in  Table  C-9  and  Fig.  C-4.  This  table  demonstrates  the 

shaping  is  not  difficult  once  the  pattern  develops  and  experience  is 
gained . 

For  example,  starting  with  the  nominal  plant  (poles  at  -0.3343  and 
-1.0054)  the  slope  of  the  magnitude  line  after  the  second  pole,  on  the 
Bode  Plot,  is  -12  db/octave.  At  this  slope  the  value  of  is  11.0  db 
at  to  ■  2  rad/sec  and  -1.0  db  at  to  =  4  rad/sec.  This  is  slightly  low  for 
the  2  rad/sec  design  bound  (approximately  13  db)  and  is  also  low  for  the 
4  rad/sec  design  bound  (approximately  2  db) .  Thus  the  slope  is 
increased  by  adding  a  zero  which  increases  the  slope  of  the  line  to  -6 
db/octave  after  the  zero  is  added.  The  problem  is  where  to  add  the 
zero.  Using  a  straight  edge  at  a  slope  of  -b  db/octave,  line  up  the  13 
db  point  at  2  rad/sec  and  the  current  magnitude  line  on  the  Bode  Plot 
whose  slope  is  -12  db/octave.  The  point  where  the  straight  edge  and  the 
-12  db/octave  line  cross,  is  the  value  oi  the  zero  which  satisfies  the 
to  «  2  rad/sec  design  bound  of  13  db.  Remember  this  is  a  straight  line 
approximation  and  the  corner  points  change  depending  on  other  roots  near 


Looking  ahead,  the  slope  of  the  line  needs  decreasing  to  achieve  a 
magnitude  close  to  the  4  rad/sec  design  bound  of  2.2  db.  Thus,  need  to 
add  another  pole  and  this  lovers  the  curve  at  the  2  rad/sec  point. 
Therefore  Instead  of  using  the  point  that  the  straight  edge  and  the 
-12  db  line  crossed  (1.7),  move  the  zero  to  1.3.  After  passing  the  2 
rad/sec  point  the  curve  needs  to  decrease  quickly  to  achieve  the  4 
rad/sec  bound,  2  db.  Before  making  more  changes  based  on  the  magnitude 
of  check  the  phase  of  the  current  design  to  the  design  bound  on  the 

Nichols  Chart  (Fig.  5-5).  The  phase  of  the  current  design  at  ca  ■  4 
rad/sec  is  -150  degrees.  On  the  Nichols  Chart  find  the  point  on  the  4 
rad/sec  design  bound  at  -150  degrees.  The  magnitude  is  7  db.  Thus  the 
bound  used  on  the  Bode  Plot  approximation  is  7  db  for  ca  *  4  rad/sec  (not 
2  db).  Using  the  same  technique  as  above,  add  a  pole  to  decrease  the 
slope,  slide  the  straight  edge  at  a  slope  of  -12  db /octave  along  the 
7  db  at  4  rad/sec  point  until  it  crosses  the  previous  -6  db/octave 
slope.  The  crossing  point  is  about  3.2,  thus  add  a  pole  at  3.2. 

The  curve  now  looks  good  until  20  rad/sec  where  the  slope  needs 
decreasing  again.  Using  the  same  technique  as  above,  add  a  zero  at 
about  12.5.  Now  check  the  current  results  on  the  Nichols  Chart  for  the 
phase  angle  and  magnitudes.  From  Fig.  5-5  the  value  for  the  current 
design  at  2  rad/sec  is  13  db  and  -139  degrees.  From  the  Nichols  Chart 
the  design  bound  magnitude  for  2  rad/sec  at  -139  degrees  is  14  db.  Thus 
the  current  design  is  about  1  db  too  low  and  requires  raising.  Notice 
that  as  the  phase  angle  becomes  more  positive  (L^q  curve  moves  left  on 
the  Nichols  Chart)  the  magnitude  of  the  design  bound  decreases,  for 
example  the  magnitude  for  the  0)  =  2  rad/sec  design  bound  at  -130  degrees 
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is  13  db.  The  value  of  the  current  design  for  w  =  4  rad/sec  is  7  db  and 
-142  degrees.  From  the  Nichols  Chart  notice  the  design  bound  is  6  db  at 
-142  degrees.  Thus  the  current  design  is  barely  satisfactory,  but  it 
violates  the  2  rad/sec  design  bound. 

Try  replacing  the  last  pole-zero  combination.  Remember  it  is 
helpful  to  add  phase  as  well  as  increase  the  magnitude.  This  is  accom¬ 
plished  by  adding  the  pole  later  (use  3.5  instead  of  3.2)  which  increases 
the  magnitude  and  adding  the  zero  earlier  (use  6.5  instead  of  12.5) 
which  adds  phase  angle.  From  the  Nichols  Chart  (Fig.  5-5  and  5-6), 
these  points  satisfy  the  bounds.  The  next  frequency  of  interest  is 
8  rad/sec.  Notice  on  the  Nichols  Chart  that  the  balance  of  frequencies 
share  the  same  design  bound.  Now  desire  the  design  to  follow  the 
8  rad/sec  bound.  Try  adding  a  complex  pole  pair  to  bend  the  curve 
around  the  design  bound  as  close  as  possible,  this  keeps  the  bandwidth 
of  the  design  as  small  as  possible.  Using  a  phase  response  chart  (Fig. 
C-7),  approximate  the  complex  pair. 

From  the  Nichols  Chart  at  10  rad/sec,  the  current  design  has  12 
degrees  of  phase  difference  between  its  value  and  the  design  bound. 

Using  the  phase  response  chart,  use  a  12  degree  angle  change  and  choose 
a  damping  ratio  of  0.7  at  10  rad/sec,  yields  a  ratio  of  frequency  to 
natural  frequency  of  0.18.  Thus  the  natural  frequency  is  56.  This 
yields  complex' poles  at  -39  +  j40  and  -39  -  j40.  Calculate  the 
values  with  these  poles  added  and  compare  it  with  the  design  bounds. 

This  now  satisfies  the  bound  (Fig.  5-6).  This  finishes  the  shaping  of 
Ljq.  Putting  together  the  winning  poles  and  zeros  yields  the 


synthesized  L  : 


20000  (s  +  1.3)(s  +  6.5) 


10  (s  -  0.3343)  (s  +  1.0054)  (s  +  39  -  j40)(s  +  39  +  j40)(s  +  3.5) 


Fig.  5-6  presents  the  open  loop  transfer  function  shaped  for  L^. 

Table  5-4  contains  the  values  at  the  frequencies  of  interest. 
Compare  these  values  with  the  Nichols  Chart  (both  phase  and  magnitude) 
as  a  check  on  the  designed.  L^q  is  in  Fig.  C-6  and  Table  C-10.  The 
does  not  meet  the  bounds  for  L^q  at  to  >  40  rad/sec  (see  Fig.  C-6). 
This  is  not  significant  because  the  lower  frequencies  are  more  important 
(designer's  discretion).  Note  that  at  this  time  in  the  design  phase, 
without  any  simulation,  the  technique  reveals  that  unless  a  change  is 
made  in  the  L^q  design,  the  system  responses  may  not  be  acceptable  at 
the  higher  frequencies  (<o  >  40) . 


Table  5-4  Values 


.  .  _ 20000  (s  +  1.3)  (s  +  6.5) _ 

10  L1  3  (s-0.3343) (s+1.0054) (s+3.5) (s+39-j40) (s+39+j40) 


Freq. 

Magnitude 

Phase 

(rad/sec) 

(db) 

(degrees) 

0.1000 

32.8692 

-165.5262 

0.2000 

31.8553 

-153.4085 

0.4000 

29.1167 

-138.0508 

0.8000 

24.0916 

-126.5848 

1.0000 

22.0762 

-124.3932 

2:0000 

15.0893 

-121.3318 

4.0000 

7.3771 

-121.6171 

8.0000 

-0.4459 

-121.4610 

10.0000 

-2.8098 

-121.7914 

20.0000 

-9.6176 

-129.7032 

40.0000 

-16.7481 

-159.1405 

80.0000 

-28.9832 

-210.3114 

100.0000 

-34.3069 

-223.4846 
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9.  The  Compensator  G 


With  L^q  formed,  it  is  easy  to  calculate  the  compensator,  since 
■  Gj  P^q.  The  open  loop  transfer  functions  for  all  three  flight 
conditions  are  now  calculated.  The  open  loop  transfer  function  for  the 
three  flight  conditions  of  the  first  loop  is  defined 

where  i  *  1,  2  or  3.  This  generates  the  open  loop  transfer  function  for 
each  flight  condition,  from  which  the  closed  loop  system  is  found. 

Tables  5-4  thru  5-7  show  the  compensator  and  open  loop  transfer 
functions. 

The  method  to  perform  the  design  is  similar  to  the  loop.  The 
results  are  in  Appendix  C. 


10.  Summary 


This  completes  the  first  part  of  the  design  which  presents  the  two 


loops,  loop  1  (L^)  and  loop  3  (L^).  The  next  chapter  presents  the 
design  of  the  last  loop,  L^.  The  last  two  chapters  present  and  discuss 


the  simulations  of  the  system.  Also,  the  design  is  compared  with  that 

* 

of  the  previous  thesis  (2)  which  uses  a  different  technique. 
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Table  5-5  Gl  Values 


Cl 

72968.9364(s  +  1.3)(s 

+  6.5) 

G1 

(s+3.5) (s+39+j40) (s+39-j40) 

Freq. 

Magnitude 

Phase 

(rad/sec) 

(db) 

(degrees) 

0.1000 

35.0559 

3.5003 

0.2000 

35.1243 

6.9517 

0.4000 

35.3857 

13.5316 

0.8000 

36.2716 

24.6034 

1.0000 

36.8123 

28.9373 

2.0000 

39.4710 

41.4688 

4.0000 

42.9981 

49.0513 

8.0000 

46.9956 

53.7688 

10.0000 

48.4810 

54.3821 

20.0000 

53.6781 

48.3766 

40.0000 

58.5797 

19.8985 

80.0000 

58.3835 

-30.7920 

100.0000 

56.9359 

-43.8691 

Table  5-6  (Lj)^  Values 

(L  .  .  _ 14887.6143  (s  +  1.3)(s  +  6.5) _ 

1  1  ”  (s+0.0811) (s+2. 1268) (s+3.5) (s+39+j40) (s+39-j40) 


Freq. 

.Magnitude 

Phase 

(rad/sec) 

(db) 

(degrees) 

0.1000 

32.4903 

-50.1496 

0.2000 

28.0434 

-66.3479 

0.4000 

22.6577 

-75.6586 

0.8000 

17.2298 

-80.2219 

1.0000 

15.5555 

-81.6086 

2.Q000 

10.3308 

-89.4492 

4.0000 

4.0262 

-101.7876 

8.0000 

-3.2314 

-110.7626 

10.0000 

-5.5178 

-113.1465 

20.0000 

-12.2184 

-125.3210 

40.0000 

-19.3214 

-156.9418 

80.0000 

-31.5496 

-209.2111 

100.0000 

-36.8724 

-222.6042 
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Table  5-7 

(L^)^  Values 

10499.2715 

(s  +  1. 3) (s  +  6.5 

(s-0. 3328) (s+1.2021) (s+3.5) (s+39+j40) (s+39-J40) 


Freq. 

Magnitude 

Phase 

(rad/sec) 

(db) 

(degrees) 

0.1000 

25.7683 

-164.5305 

0.2000 

24.7845 

-151.4902 

0.4000 

22.1654 

-134.6338 

0.8000 

17.4858 

-121.6278 

1.0000 

15.6336 

-119.2265 

2.0000 

9.1321 

-116.9707 

4.0000 

1.6705 

-118.9780 

8.0000 

-6.0722 

-120.0678 

10.0000 

-8.4258 

-120.6694 

20.0000 

-15.2197 

-129.1371 

40.0000 

-22.3467 

-158.8568 

80.0000 

-34.5809 

-210.1695 

100.0000 

-39.9045 

-223.3710 

11  22  33  q12  q23  q31  q13  q32  Q21 


This  loop  contains  less  overdesign  than  the  first  two  designs 
because  all  terms,  except  the  L2  terms,  are  known  exactly.  Now  define 
an  effective  plant  (the  plant  the  system  sees)  Q22e  as: 


_ ^22 _ 

r12  (1  +  L3)  +  Y23  (1  +  Lx) 
(1  +  L^d  +  L3)  -  y13 


(6-2) 


or  rewriting  in  an  equivalent  form: 

q00[(l  +  L, ) ( 1  +  L_) 


22e  [(1  +  L^d  +  L3)  -  y13J  "  LY12  U  +  L3)  +  Y23  (1  +  Lj)  -I*] 

* 

Substitute  Q22e  into  Eq.  6-1  and  rewrite  as: 


F  L 
=  22  20 

22  (1  +  L2Q) 


where  L2Q  -  G2g  Q22e  «  G2  Q22»  Use  Q22g  to  design  L2Q  and  then  derive 
the  compensator,  G2>  from  L2q. 

This  has  the  same  form  as  the  S1S0  systems  of  Chapter  11. 

Generate  the  plant  templates  using  Q  as  the  plant  (as 


values 


done  for  L.  and  L.  in  Section  V-6).  Table  6-1  contains  the 

1  3  22e 

and  the  templates  are  in  Fig.  6-1.  Appendix  D  contains  the  Q£2e  calcula¬ 
tions. 


Table  6-1  Q„„  Values 

22e 


Freq. 

FC 

//I 

FC 

n 

FC 

it  3 

(rad/ sec) 

(db) 

(degrees) 

(db) 

(degrees) 

(db) 

(degrees) 

0.1 

14.65 

11.98 

12.72 

21.36 

24.48 

87.27 

0.2 

15.90 

17.09 

14.62 

19.95 

26.77 

118.18 

0.4 

17.73 

18.48 

16.02 

17.45 

27.32 

-220.24 

0.8 

20.49 

12.92 

18.58 

15.97 

25.96 

-203.39 

1.0 

21.37 

5.94 

19.99 

11.33 

25.51 

-196.86 

2.0 

22.93 

-17.96 

25.61 

-23.91 

24.33 

-183.65 

4.0 

24.93 

-110.00 

15.63 

-157.93 

19.94 

-180.43 

8.0 

7.72 

-158.38 

0.076 

-168.43 

11.39 

-178.00 

10.0 

3.27 

-160.53 

-4.13 

-168.80 

8.08 

-177.15 

20.0 

-9.27 

-159.08 

-16.48 

-165.83 

-3.05 

-174.62 

40.0 

-20.64 

-148.40 

-28.13 

-156.39 

-14.75 

-170.14 

80.0 

-30.42 

-130.91 

-38.63 

-139.97 

-26.38 

-161.38 

100.0 

-33.14 

-124.94 

-41.62 

-133.75 

-30.02 

-157.26 

This  loop  requires  an  upper  and  lower  bound,  generated  by  using  b^ 
and  a^2  (Table  5-2  and  Fig.  5-3).  Calculate  the  difference  between  the 
two  bounds  at  the  frequencies  of  interest,  call  this  difference  A  (Table 
6-2).  The  design  proceeds  as  explained  in  Section  11-9  with  A  defined 
as  the  maximum  change  allowed  at  each  frequency.  Slide  the  templates 
along  the  Nichols  Chart  to  achieve  these  maximum  changes.  Transfer  the 
reference  point  from  the  template  to  the  Nichols  Chart  at  each  point 
where  the  template  exactly  satisfies  A  for  that  phase  angle  (as  done  in 
Section  11-9).  Do  this  at  enough  phase  angles  to  enable  drawing  a 
representative  curve  at  each  frequency.  Remember  the  UHF  bound  may  not 


Table 

6-2  (A  Values, 

A  -  |b22|  - 

la22l) 

4 

4 

« 

Freq. 

b22 

a22 

A 

(rad/sec) 

(db) 

(db) 

(db) 

« 

0.1 

0.0063 

-0.036 

0.043 

0.2 

0.13 

-0.14 

0.27 

. 

0.4 

0.62 

-0.54 

1.16 

0.8 

2.35 

-1.98 

4.33 

1.0 

3.25 

-2.91 

6.16 

2.0 

-1.10 

-8. 15 

7.05 

4.0 

-9.81 

-17.15 

7.34 

* 

8.0 

-15.69 

-28.68 

12.99 

10.0 

-17.38 

-32.78 

15.40 

20.0 

-22.81 

-46.86 

24.05 

40.0 

-28.61 

-63.02 

34.41 

80.0 

-34.56 

-80.44 

45.88 

» 

100.0 

-36.49 

-86.18 

49.69 

* 

be  violated.  For  several  of  the  templates  (<a  >.  8  rad/sec)  the  UHF  Bound 
Is  the  nominal  bound.  Other  templates  (ta  <  8  rad /sec)  use  the  A  limit 
as  the  nominal  bound  except  at  the  phase  angles  where  the  UHF  Bound  is 
more  restrictive.  After  using  all  of  the  templates,  the  design 
bounds  are  complete  (see  Fig.  6-2). 

The  method  used  to  obtain  the  bounds  for  co  £  0.4  rad/sec  is  an 

analytical  method.  The  templates  are  very  large  and  the  differences 

allowed  very  small.  The  template  graphically  finds  a  point,  call  this 

L  ,  such  that  the  difference  between  L  and  a  point  on  the  template,  and 

and  any  other  point  on  the  template  is  less  than  or  equal  to  A.  Thus 

take  two  points  on  the  template  and  find  a  point  at  each  frequency  of 

interest  such  that  the  difference  between  the  distance  from  L  to  one 

a 
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point  of  the  template  and  the  distance  from  and  the  other  point  is 

less  than  or  equal  to  A .  In  this  thesis  this  is  relatively  simple  (see 
Appendix  D) .  Observe  the  results  from  the  other  loops,  the  pattern 
indicates  the  bounds  for  £.0.8  rad/sec  is  greater  than  30  db  (Fig. 
6-2).  The  bounds  are:  48  db  at  gj  =  0.1  rad/sec,  32  db  at  gj  ■  0.2 
rad/sec  and  33  db  at  u  ■  0.4  rad/sec  (see  Appendix  D  for  the 
calculations) . 

A  design  bound  now  exists  at  all  the  frequencies  of  interest. 

Shape  The  procedure  is  the  same  as  the  shaping  performed  for  the 

other  two  loops  in  Chapter  V  (see  also  Section  11-11).  The  results  are 
in  Fig.  6-3  and  Table  6-3. 


Table  6-3  L„„  Values 


(s+0.9848) (s+0. 1) (s+0. 1) (s+O.4) (s+0. 8) (s+8) (s+10) 

(s+636-j848) (s+636-j848) 


Freq. 

(rad/sec) 

(db)  120 

(degrees) 

0.1 

47.20 

-69.01 

0.2 

40.93 

-93.47 

0.4 

32.81 

-107.27 

0.8 

23.92 

-105.68 

1.0 

21.27 

-101.26 

2.0 

15.11 

-82.72 

4.0 

11.52 

-75.22 

8.0 

7.45 

-86.91 

10.0 

5.58 

-92.75 

20.0 

-1.90 

-106.38 

40.0 

-9.95 

-107.69 

80.0 

-17.01 

-104.71 

100.0 

-19.10 

-104.36 

V 


With  L 2q  designed,  find  (L^  using  L2  *  G2  ^see  Fig*  • 

The  prefilter  simply  moves  the  open  loop  transfer  function,  L2>  into  the 
acceptable  region  defined  by  the  upper  and  lower  bounds  (Fig.  5-3).  The 
L 2Q  design  only  guaranteed  the  system  responses  do  not  exceed  the 
maximum  change,  A.  Compare  the  results  in  Table  6-4,  the  closed  loop 
system  [I^/O  +  I^)],  with  the  values  of  b^  and  a22*  The  va^ues  °f  the 
closed  loop  system  is  within  the  acceptable  region  for  0)  £  1.0  rad/sec. 
For  the  frequencies  above  1.0  rad/sec  the  values  need  decreasing. 


Table  6-4 

(  L2  ) 

Closed  Loop  System 

/Jjl) 

Values  (x  +  j) 

/J2L\ 

Maximum 

l  i+s/i 

l 1+L2  k 

Difference 

Freq. 

(rad/sec) 

(db) 

(db) 

(db) 

(db) 

0.1 

-0.014 

-0.015 

-0.00018 

-0.015 

0.2 

0.004 

0.030 

0.020 

0.026 

0.4 

0.06 

0.20 

0.11 

0.14 

0.8 

0.13 

0.68 

0.43 

0.55 

1.0 

0.12 

0.79 

0.57 

0.67 

2.0 

-0.32 

-0.59 

0.23 

0.82 

4.0 

-0.81 

-2.62 

-0.80 

1.82 

8.0 

-0.88 

-4.16 

-1.28 

3.28 

10.0 

-0.89 

-5.06 

-1.51 

4.17 

20.0 

-2.67 

-11.10 

-5.14 

8.43 

40.0 

-9.53 

-19.33 

-13.02 

9.80 

80.0 

-16.88 

-26.64 

-20.32 

9.76 

100.0 

-18.91 

-28.65 

-22.45 

9.74 

Design  the  prefilter  f22*  to  shift  the  response  magnitude  into  the 
acceptable  region.  Use  only  the  magnitude  values  because  the  design 


assured  the  appropriate  phase  angle.  First  draw  the  envelope  into  which 


the  prefilter  needs  to  exist  (Fig.  6-4).  Do  this  by  comparing  the 


largest  value  of  the  closed  loop  system,  ,  to  the  upper  bound,  T^, 

and  the  lowest  value  of  the  closed  loop  system,  T  to  the  lower 
bound,  T^,  at  all  frequencies  of  interest  (Fig.  2-9).  These 
differences,  which  must  be  greater  than  or  equal  to  zero,  define  a  range 
on  the  Bode  Plot  in  which  the  prefilter  exists  such  that  b^  and  a ^ 
envelope  the  system  response  curves  (Fig.  6-4).  For  example,  at  u  »  2 
rad/sec,  the  largest  value  from  the  closed  loop  system  responses  is  0.23 
db  (FC  #3)  and  the  lowest  value  is  -0.32  db  (FC  #1)  (Table  6-4).  The 
upper  bound,  is  -1.10  db  and  the  lower  bound,  a^2  is  -8.15  db  (Fig. 
5-3  and  Table  6-2).  Therefore  the  upper  limit  of  F 22  is  -1.10  -  0.23  * 
-1.33  db  and  the  lower  limit  is  -8.15  -  -0.32  ■  -7.83  db.  Notice  a  line 
sloped  at  approximately  6  db/octave  satisfies  the  bounds  between  1  and 
100  rad/sec.  Therefore  try  a  simple  prefilter  with  one  pole.  Remember, 
the  prefilter  values  are  very  small  for  0)^1  rad/sec  because  the  A  at 
those  frequencies  are  very  small.  The  prefilter  chosen  (see  Table  6-5) 
is: 


F _ hm. 

22  (s  +  1.4) 


This  concludes  the  design.  The  four  parts  of  the  system,  the  three 
compensators  (one  for  each  loop)  and  the  prefilter  for  loop  2,  designed 
are: 


=  72968.936  (s  +  1.3)(s  +  6.5) 

°1  (s  +  3 . 5) (s  +  39  +  j 40)  (s  +  39  -  j40) 
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G  =  2521167.3  (s  +  0.18)(s  +  l)(s  +  1) (s  +  1.5)(s  +  1.2)(s  +  30) 

1  (s+O. 1) (s+0.1) (s+0.4) (s+0.8) (s+8) (s+10) (s+636+j848) (s+636-j848) 


G  „  _ 98.917 _ 

3  s(s  +  0.7  +  j  0. 7 14) (s  +  0.7  -  j  0.714) 


F  -  _ 

22  (s  +  1.4) 


Table  6-5  F ^  Values 

r  ,  1‘404 

*22  (s  +  1.4) 

Freq. 

(rad/sec) 

^22 

<db) 

(degrees) 

0.1 

0.00 

-4.09 

0.2 

-0.07 

-8.13 

0.4 

-0.32 

-15.95 

0.8 

-1.20 

-29.74 

1.0 

-1.77 

-35.54 

2.0 

-4.81 

-55.01 

4.0 

-9.60 

-70.71 

8.0 

-15.25 

-80.07 

10.0 

-17.14 

-82.03 

20.0 

-23.10 

-86.00 

40.0 

-29.10 

-88.00 

80.0 

-35.12 

-89.00 

100.0 

-37.06 

-89.20 

4.  Summary 


This  completes  the  lateral  controller  designed  Lor  this  thesis. 
The  last  two  chapters  walked  through  the  design  technique  used  in  this 
thesis.  The  next  few  chapters  present  and  discuss  the  simulations  of 


VII.  SIMULATION  RESULTS 


1.  Introduction 

This  chapter  presents  the  responses  of  the  lateral  controller 
designed  in  Chapters  V  and  VI.  First  the  three  compensator  designs, 

G^,  and  G^  are  presented.  Then  the  three  responses,  t^»  ^2*  ant*  C32 
at  all  three  flight  conditions  are  presented,  followed  by  the  three 
responses  at  different  magnitude  step  inputs. 

2.  Compensator  Characteristics 

From  the  original  specifications  in  Chapter  V  both  loop  1  and  loop 

3  response  steady  state  values  are  zero.  Thus  the  denominator  of  the 

compensators  must  contain  a  root  (pole)  at  zero.  The  loop  3  transfer 

function  contains  a  pole  at  zero,  but  the  loop  1  transfer  function  does 

not.  Therefore  a  pole  at  zero  is  added  to  the  loop  1  transfer  function 

Also  a  root  is  added  to  the  numerator  (called  a  zero)  such  that  the 

response  characteristics  remain  relatively  unchanged  at  the  frequencies 

of  interest.  The  modified  G.  is  called  G, 

1  lmod 

G  =  72968.9364  (s  ♦  0.25)(s  +  1.3)(s  +  6.5) 
lmod  s (s  +  3.5)  (s  +  39  +  j40)(s  +  39  -  j40) 

The  two  responses  are  presented  in  Fig.  7-1. 

Basically,  bandwidth  indicates  to  which  frequency  the  transfer 
function  amplifies  its  input  signals.  The  bandwidth  of  the  system  is 
the  range  of  frequencies  at  which  the  magnitude  of  the  function  is 
greater  than  0  db.  Input  signal  attenuation  is  desired  as  soon  as 


possible . 


The  bandwidth  of  the  compensators  are  important  because  it  tells 

for  which  frequencies  the  compensator  inputs  are  amplified  (in  real 

systems  this  input  signal  might  be  corrupted  by  noise  and  other  unwanted 

signals).  The  responses  of  the  four  compensators  G,,  G,  G„  and  G„ 

1  imod  2  J 

are  in  Figs.  7-1.  Remember  these  values  are  in  radians  and  rad/sec. 

Using  the  units  of  degrees  and  hertz,  the  bandwidths  of  G^  and  G^  are 

reasonable,  though  the  bandwidth  of  G^  and  G^mod  are  large.  The 

bandwidth  of  G,  ,  is  72000  rad/sec  (1250  hz).  G ,'s  bandwidth  is 
lmod  1 

between  0.1  and  1250  hz.  The  bandwidth  or  G^  is  1610  rad/sec  (0.08  hz) . 
The  bandwidth  of  G^  is  4.6  rad/sec  (1.1  hz). 

Notice  the  bandwidth  of  the  compensator  with  the  commanded  input, 
G^,  is  much  lower  than  even  G^,  though  G^  is  designed  last. 

3.  Loop  Responses 

To  simulate  the  loop  responses,  use  the  system  equations. 

Reviewing  the  original  3x3  M1M0  system,  the  matrix  equation  is; 

T  =  [1  +  GP]"1  PFGR  (7-1) 

There  are  nine  terms  in  "T.  Only  three  of  the  terms  are  nonzero  because 
only  loop  2  has  an  input  command.  Thus  the  six  terms  which  have  an 
input  command  to  loop  1  or  3  ( t £ ^ *  t31’  C13’  C23^  are  zero •  Solve  the 
matrix  equation  for  the  three  scalar  equations  of  interest.  The 


equations  are; 


Solve  Eq.  7-2b  in  terms  of  and  (see  Eq.  6-1).  Then  solve  Eqs. 

7-2a  and  c  in  terms  of  t^2  (see  Appendix  E) . 

Generate  the  loop  response  using  the  loop  equations  derived.  Fig. 

7-2  through  7-5  contain  the  design  response  bounds.  Figs.  7-6  contains 

the  FC  #2  responses  to  an  input  of  1  ft/sec.  The  time  specifications 

for  all  of  the  flight  conditions  and  the  design  bounds  are  in  Table  7-1 

Comparing  the  response  specifications  (Table  7-1),  all  of  the 

responses  satisfy  the  design  bounds.  Note  that  t^  response  magnitudes 

are  about  an  order  of  magnitude  less  than  the  design  bounds.  Thus,  the 

tp  being  larger  than  the  design  bounds  t  (remember  in  Section  V-4,  t 

is  not  one  of  the  important  specification  criteria)  is  not  a  problem. 

Note  also  the  shapes  of  the  curves  for  the  bounds  and  responses  are 

-14 

similar  (notice  the  magnitude  of  the  Fig.  7-7 (a)  curve  is  10 
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Table  7-1  Loop  Responses 


(.Step 

input  of  1 

ft/sec) 

t 

t  , 

t 

t 

M 

F 

r 

d 

P 

s 

P 

V 

b12 

1.60 

0.0499 

0.00 

(ti2)i 

6.22 

0.0034 

0.00 

“llh 

15.97 

0.000 

0.00 

2.930 

14.87 

0.033 

0.00 

b22 

0.751 

0.950 

1.86 

5.58 

1.31 

1.00 

(t22)i 

0.89 

11.55 

16.17 

1.56 

1.01 

1.00 

<t22>2 

1.46 

2.11 

3.08 

4.61 

1.09 

1.00 

<*«>, 

0.98 

1.37 

2.51 

1.22 

1.00 

a22 

2.12 

3.76 

1.00 

1.00 

b32 

0.237 

0.309 

0.00 

(t32)i 

0.30 

0.0034 

0.00 

(t32)2 

2.40 

0.0032 

0.00 

(t32J3 

0.810 

0.000935 

0.00 

4.  Comparison  of  Loop  Responses  and  Bounds 

The  comparison  of  specifications  is  done  in  the  time  domain  since 
the  specifications  are  defined  in  the  time  domain.  Table  7-1  lists  the 


time  response  specifications  for  the  t^»  ^2  an<*  C32  resPonses  an<* 
their  respective  bounds  (b^*  ^22*  anc*  ^32^'  ^emember  the  bounds  are 
derived  using  only  a  few  of  the  criteria  as  important  (t^) ,  ^  and  F^, 


Section  V-4) . 


All  of  the  magnitude  responses  are  below  the  design  bound,  b^- 
Mp  for  b^  is  0.05  (.step  input  =  1  it/sec)  and  all  three  of  the  flight 
conditions  are  below  this.  The  t ^  responses  are  within  their  defined 
envelope  (all  three  flight  conditions  have  a  tf  value  between  the  b ^ 
and  a„„  values  of  0.751  and  2.12,  as  well  as  a  M  between  1.0  and  1.3; 

22  p 

Table  7-1).  The  t.„  responses  satisfy  the  b.»  bound  (b._  bound  M  ■  0.3 
32  J2  32  p 

and  the  three  flight  conditions  are  almost  two  magnitudes  smaller). 

All  three  loop  transfer  functions  designed  satisfy  the  design 
specifications  for  all  of  the  flight  conditions.  Thus  the  design  is 
acceptable. 


5.  Different  Command  Inputs 

Now  try  several  different  magnitude  inputs.  The  inputs  used,  or 
the  side  velocity  commanded,  are  22.4  ft/sec,  32  ft/sec  and  48  ft/sec. 
Only  a  representative  response  (FC  #2)  for  each  input  is  shown  here. 
The  results  are  in  Figs.  7-7  through  7-9  and  Appendix  E. 


All  of  the  responses  retain  their  original  t  ,  t  ,  t ,  and  t  .  The 

p  s  d  r 


F^  changes  appropriately  for  loop  2,  to  match  the  commanded  input,  i.e. 


for  a  step  input  of  32  ft/sec  the  F^  is  32  ft/sec  (Table  7-3  and  Fig. 


7-6b  and  Appendix  E) .  The  F^  of  the  other  two  loops  are  zero.  The  M^, 


basically  a  measurement  of  the  amount  of  overshoot,  remained  linear.  M 


remained  at  a  constant  percentage  overshoot  in  loop  2.  For  example,  for 


an  input  of  1  ft/sec,  =  1.09  ft/sec  and  for  an  input  of  32  ft/sec. 


=  34.86  ft/sec  (FC  it  2) .  Both  responses  have  an  overshoot  of  9%.  An 


input  of  1  ft/sec  to  loop  1  results  in  »  0.0034  rad  (FC  01)  and  an 


input  of  32  ft/sec  results  in  M  =  0.108  rad  (32  x  0.0034).  An  input  of 


r(t)  -  22.4  ft /sec 
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Table  7-3  Loop  2  Responses  to  Different  Inputs 


v 

V 

j 

'j 

v 

\ 


t 

t  , 

t 

t 

M 

r 

d 

P 

s 

P 

Step  Input 

of  22.4 

f t/sec 

0.751 

0.950 

1.86 

5.58 

29. 

0.87 

11.55 

16.17 

1.56 

22. 

1.46 

2.11 

3.08 

4.61 

24. 

0.98 

1.37 

2.51 

very  large 

27. 

2.12 

3.76 

22. 

Step 

Input  of  32.0 

ft/sec 

0.751 

0.950  1.86 

5.58 

41. 

0.87 

11.55 

16.17 

1.56 

32. 

1.45 

2.11 

3.08 

4.61 

34. 

0.98 

1.37 

2.51 

very  large 

39. 

2.12 

3.76 

32. 

Step  Input  of  48.0  ft/sec 


b22 

0.751 

0.950 

1.86 

5.58 

62. 

(t22>i 

0.87 

11.55 

16.17 

1.56 

48. 

u22>2 

1.46 

2.11 

3.08 

4.61 

52. 

(t_„) _ 

0.98 

1.37 

2.51 

very  large 

59. 

Table  7-4  Loop  3  Responses  to  Different  Inputs 


cd 

t 

P 

t 

s 

M 

P 

Step 

Input  of  22.4 

f t/sec 

0.237 

very  large 

0.692 

0.30 

very  large 

0.12 

2.40 

very  large 

0.0710 

0.810 

very  large 

0.0209 

Step  Input  of  32.0 

ft/sec 

0.237 

very  large 

0.988 

0.30 

very  large 

0.17 

2.40 

very  large 

0.101 

0.810 

very  large 

0.0299 

Step  Input  of  48.0 

ft/sec 

0.237 

very  large 

1.48 

0.30 

very  large 

0.25 

2.40 

very  large 

0.152 

Table  7-5  Comparison  of  Responses 
Previous  Thesis  and  This  Thesis  (2:62-63) 


fc  n 

Step  input  of 
48  ft/sec 

FC  #2 

Step  input  of 
32  ft/sec 

FC  #3 

Step  input  of 
22.4  ft/sec 


FC  #1 

Step  input  of 
48  ft/sec 

FC  t2 

Step  input  of 
32  ft/sec 

FC  ii  3 

Step  input  of 
22.4  ft/sec 


t 

r 


Previous  Thesis 


t 

s 


1.5  sec 


7.0  sec 


0.9  sec 


5.8  sec 


0.8  sec 


3.0  sec 


This  Thesis 

0.87  sec  1.6  sec 


1.45  sec 


4.6  sec 


0.98  sec 


the  third  flight  condition  but  the  specifications  are  satisfied  by  the 


responses. 

6.  Summary 

The  compensator  responses,  G^,  and  G^  are  presented.  The 

loop  responses  are  presented  with  a  comparison  to  their  design  bounds. 
Then  a  comparison  is  made  between  the  previous  thesis's  results  (2)  and 
the  results  from  the  lateral  controller  designed  in  this  thesis. 


VIII.  DISCUSSION  AND  CONCLUSIONS 


1.  Summary  of  Thesis 

This  thesis  reviewed  the  Quantitative  Feedback  Theory  (QFT)  techni¬ 
que  used  to  design  a  lateral  controller  for  a  hypothetical  fighter  type 
aircraft.  Chapter  II  and  Appendix  A  explains  the  QFT  technique  for  a 
single  input-single  output  (S1S0)  system.  Chapter  III  explains  the 
technique  for  a  multiple  input-multiple  output  (MIMO)  system.  The 
three-by-three  aircraft  model  is  reviewed  in  Chapter  IV  and  Appendix  B. 
The  design  for  the  first  two  of  the  three  loops  designed  are  in  Chapter 
V  and  Appendix  C.  The  final  loop  design  is  in  Chapter  VI  and  Appendix 
D.  Chapter  Vll  and  Appendix  E  contain  the  system  responses  generated 
from  the  computer  simulation. 

2.  Discussion  of  the  Results 

From  the  results  in  Chapter  VII  and  Appendix  E,  the  designed 
lateral  controller  satisfies  all  of  the  design  specifications.  The 
designed  system  responses  are  better  than  the  design  specification. 

Comparing  the  two  methods,  the  previous  thesis  results  and  this 
one,  the  responses  are  similar.  This  was  expected  because  the  response 
of  Che  previous  thesis  is  used  to  generate  the  design  specifications  for 
this  thesis.  The  two  methods  do  differ.  The  Porter  technique  must  test 
each  design  for  satisfaction  of  the  performance  specifications  bounds 
and  for  robustness.  The  QFT  technique  guarantees  that  any  uncertainty 
designed  for  is  satisfied  by  the  design  solution,  thereby  the  perfor¬ 
mance  specifications  are  satisfied  and  the  robustness  (uncertainty 
range)  is  satisfied.  Remember  the  uncertainty  and  performance 


specifications  are  used  in  the  QFT  technique  to  synthesize  the  controller 
Porter's  method  used  straight  gains  in  the  feedback  loop  to  achieve  a 
compensator  to  satisfy  the  specifications.  This  generates  a  compensator 
bandwidth  of  unknown  size.  Since  the  QFT  technique  is  performed  in  the 
frequency  domain,  the  bandwidth  can  be  controlled.  Porter's  technique 
might  find  a  solution  faster  than  the  QFT  technique  if  a  solution  is 
found  in  the  first  few  iterations.  If  the  specifications  are  difficult, 
or  impossible  though,  the  Porter  method  could  take  much  longer.  The  QFT 
technique's  insight  helps  here.  The  difficult  or  impossible 
specifications  are  identified  by  the  designer  who  then  modifies  the 
design  bounds  and  specifications  if  necessary. 


3.  Possible  Improvements  for  the  Design 

There  are  several  ways  to  improve  this  design.  Do  a  design  trade¬ 
off  during  the  generation  of  the  nominal  bounds,  rearrange  the  order  in 
which  the  design  is  performed,  use  the  improved  method  or  use  a  non¬ 
diagonal  G  or  prediagonalized  system. 

A  trade-off  can  be  made  during  the  generation  of  the  nominal 
bounds.  If  one  loop  has  an  unacceptable  high  bandwidth,  and  the  other 
loops  easily  meet  their  specifications,  then  by  requiring  the  loops 
which  easily  meet  their  specifications  to  meet  more  restrictive  specifi¬ 
cations,  it  might  help  reduce  the  bandwidth  of  the  loop  with  too  high  a 
bandwidth.  This  is  done  by  decreasing  the  appropriate  design  bound 
during  the  nominal  bound  calculation.  For  example  in  Eq.  5-1: 


(5-1) 
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1 1 1 1  +  b 
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12  1 


‘13' 
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If  the  Lj  bound  is  hard  to  satisfy,  but  the  bound  is  not,  then  by 
decreasing  the  b^  value  (making  it  more  restrictive)  decreases  the 
bound  on  Lj . 

If  during  the  design,  the  designer  notices  that  a  Iood  other  than 
the  one  suspected  hardest  to  satisfy  is  hardest,  rearrange  the  order 
in  which  the  loops  are  designed.  In  this  paper,  looking  at  the  resulting 
responses,  the  order  of  design  should  change  from  performing  loop  1  and 
3  first  and  then  loop  2  to  performing  loop  3  and  2  first  and  then  loop 
1.  Loop  1  has  a  very  high  bandwidth,  therefore  by  doing  this  loop 
last,  the  bandwidth  could  be  lowered  a  little  because  there  is  less 
overdesign  in  the  last  loop  designed  (Section  V-5). 

Another  improvement  is  to  use  the  improved  method  (Section  III-5). 

This  does  not  change  the  first  loop  design.  But  the  second  loop 
designed  uses  the  loop  previously  designed,  thus  less  uncertainty  exists 
in  the  terms.  This  leads  to  less  overdesign  throughout  the  designs. 

Another  improvement  is  to  use  a  nondiagonal  G  matrix.  It  involves 
trying  to  prediagonalize  the  plant.  From  the  loop  transmission 
synthesized,  a  is  derived  which  effectively  includes  the  diagonalizing 

matrix,  H,  which  is  used  to  generate  the  diagonalized  plant  matrix,  P 


!  !  1 


H  =  P_1A 
o 


P  -  PE  »  PP  A 
eff  o 


4.  Comments  on  the  QFI  Technique 

The  QFT  technique  has  several  advantages.  Probably  the  most 
important  advantage  is  its  transparency.  The  designer  knows  during  the 
design  phase  how  well  the  design  is  doing.  He  knows  if  the  design  meets 
its  goals  or  not.  If  it  does  not  meet  its  goals,  the  technique  guides 
the  designer  to  the  source  of  the  problem,  and  thus  reveals  possible 
corrective  actions. 

Another  advantage  is  the  algebra  is  relatively  simple.  There  are 
no  complex  integrals  or  derivatives  to  manage.  This  makes  the  mechanics 
of  the  technique  easy  and  the  programming  (.except  for  the  graphical 
part)  of  the  problem  simple. 

Probably  the  biggest  disadvantage  is  the  amount  of  overdesign. 

There  often  is  overdesign.  In  other  words,  though  the  design  is 
guaranteed  to  meet  the  specifications,  more  bandwidth  is  used  than  is 
necessary.  The  amount  more  than  necessary  cannot  be  determined  without 
some  analysis  of  the  design.  Reference  B  discussed  means  for  reducing 
the  overdesign. 

A  point  of  argument  is  the  amount  of  time  and  effort  needed  o 

design  a  M1M0  system  because  of  the  number  of  loops  necessary  to  design. 

2 

The  QFT  technique  reduces  the  n  x  n  MIMO  problem  into  n  equivalent  SISO 
problems,  with  n  loops  to  solve.  Consider  the  alternatives  though. 

Huge  matrices  with  many  assumptions  to  simplify  the  matrices  or  fill  out 


the  matrices,  requiring  a  computer  to  solve. 


On  the  other  hand,  the  QFT  technique  requires  the  design  of  many  S1S0 
problems.  Using  a  computer  program  (other  techniques  require  a  computer 
program  for  almost  any  real  problem),  the  SISO  solutions  are  found  very 
quickly  and  remember  the  solution  is  guaranteed  for  minimum  phase  plants 
to  meet  the  design  criteria  (the  design  could  be  also  performed 
manually).  If  the  n  x  n  M1M0  problem  does  not  have  certain  inputs  or 
coupling  terms,  the  number  of  equivalent  problems  to  solve  might  drop 
significantly.  For  example,  in  this  thesis,  only  loop  2  had  a  non-zero 
command  input.  Thus  only  three  of  the  possible  nine  equivalent  SISO 
problems  from  the  3x3  MIMO  system  needed  solving.  Many  of  the  other 
techniques  must  still  solve  huge  matrices  unless  the  matrices  are 
reconfigured  which  could  take  a  lot  of  work,  time  and  approximations. 

5.  Recommendations  and  Conclusions 

The  QFT  technique  is  a  method  which  needs  further  investigation. 

The  technique's  theory  says  it  could  work  on  most  any  minimum  phase 
problem,  linear  or  non-linear.  Only  some  of  these  problems  have  been 
investigated.  This  is  one  of  several  theses  from  AF1T  investigating  the 
simple  linear  minimum  phase  problem.  This  is  the  first  AF1T  thesis  to 
investigate  a  system  higher  than  a  two-by-two.  The  technique  does  get 
more  involved  with  each  additional  loop,  which  is  expected.  Throughout 
these  investigations  the  QFT  technique  has  been  demonstrated  for  this 
type  of  problem.  Now  non-linear  and  non-minimum  phase  systems  need 
investigating  more  fully  (some  investigation  has  been  done  already) . 
Another  area  is  to  complicate  the  plant  by  adding  bending  modes  and 
coupling  the  lateral  and  longitudinal  axes.  A  computer  program  (CAT 
system)  with  a  seieen  showing  the  bode  1’lot  and  Nichols  Chart 


vw: 


1  15 


.*  ,*  f 


together,  as  well  as  a  time  response  plot  would  be  very  helptul  and 
speed  up  the  method  significantly. 

More  research  is  needed  to  find  the  limits  of  this  technique  so  the 
designer  can  put  this  method  in  his  repertoire  and  use  it  to  its  best 
advantage. 

In  conclusion,  the  technique  works  as  advertised,  the  design  met 
and  exceeded  all  of  the  design  specifications.  The  technique  uses 
simple  algebra  and  a  graphical  technique,  which  makes  it  transparent 
from  the  solution  back  to  the  original  specifications. 
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Appendix  A 

Loop  Transmission  Shaping  Example 


This  appendix  contains  an  example  of  shaping  a  loop  transmission. 
This  is  a  copy  of  a  handout  Dr.  Horowitz  distributed  in  a  class  at  AF1T 
The  example,  which  is  complete  within  itself,  is  meant  to  give  the 
reader  a  better  understanding  of  how  to  shape  a  loop  transmission. 


Example :  Shaping  of  a  nominal  loop  transmission  L^(jai)  to  satisfy 

boundaries  B(uO  on  Nit  hoi s  Chart. 

Previous  notes  have  described  how  tolerances  on  the  closed-loop  system 

frequency  response  are  readily  translated  into  bounds  on  a  nominal  loop 

transmission  function  L^fjca).  In  Fig.  1,  for  example,  ( j  2 )  must  be  on 

or  above  the  curve  labelled  B(2),  etc.  B,  is  the  "universal  high-f requenev 

h 

boundary"  applicable,  in  this  example,  to  i.e.  L^(juj)  (for  ujJ^U) 

must  be  contained  in  the  closed  curve  B^  in  Fig.  1.  Additional  specifica¬ 
tion  is  eL'^»  where  is  excess  of  poles  over  zeros  of  L^Cs).  Also,  is 
to  be  Type  1  (one  pole  at  the  origin).  We  proceed  to  describe  a  reasonable 
procedure  for  choosing  a  rational  function  Lq(s)  which  satisfies  the  above 
specif ications. 

In  our  first  step,  we  try  to  find  the  B(w)  which  "dominates"  L^(ja')- 

E.g.  suppose  LQ1 (j4)=Odb/-135°  (point  k  in  Fig.  1).  But  at  m=l,  |  I-q ( j ) 1 

needed  is  =27db.  In  order  to  decrease  | |  from  27db  to  about  Odb  i;.  2 
2 

octaves  (4/1  =  2  ),  the  oi  J  ( joo) |  would  have  to  be,  on  the  average, 

about  -14db/octave ,  involving  ^^<-180°.  We  assume  "absolute"  stability 
is  required  here  for  Lg(jw)  with  a  margin  of  40”,  not  just  at  crossover 
("crossover"  is  defined  as  the  frequency  at  which  | L^( j u>) | -1) .  Hence  B(l) 
dominates  Lq(Juj),  at  least  more  than  B(4).  In  the  same  way  we  see  that 
B ( 1 )  dominates  over  all  other  B(uj)  in  Fig.  1. 

The  B(w)  for  uj<1  are  not  shown  in  Fig.  1.  Wc  shall  assume  that  for 
liXl,  a  slope  of  -6  db/octavo  (with  27db  at  ic=l,  i.e.  33  db  at  u  = .  5 ,  39db 


at  ia“.25  etc.),  suffices.  Wc  can  tolerate  =—140°  for  wjl ,  so  we  choose  a 
lag  corner  frequency  (symbol  lacf )  at  up  1  (i.e.  pole  at  -1),  and  set 

I  h  ( J  a')  |  (asymptotic)  at  30  db  (to  allow  for  the  -3db  n  r-v:  imo.  1;  . 
our  Lq  is  so  far:  L  = 31 . 6 /s ( s+l ) , 

Fig.  2. 

A-2 


whose  phase  J\  ( j <--)  is  shot  died  ;  a 

_ 


Lg^(jiii)  violates  tlie  -140*  bound  at  w>1.2,  so  a  lead  corner  frequent-;. 


(symbol  lccf )  is  needed.  Where  should  it  be  located?  At  ^-5,  /Lqj  16' 

(see  Fig.  2),  so  29°  lead  is  needed;  but  we  know  that  later  there  will  bo 

a  second  lac  f .  so  allow  say  additional  15°  for  it  giving  15+29~45°  lead 

required  at  w=5;  which  is  achieved  by  a  lecf  at  u)=5  i.e.  a  zero  at  -5. 

The  resulting  31 . 6  (1+  -|) ,  whose  phase  ^ l^O")  is  sketched  in  Fig.  2 

s (14s) 

In  the  Nichols  Chart’  (N.C.)  we  are  (=-<a“lO  or  so)  in  the  region  where  the 

maximum  phase  lap  allowed  is  135°  (i.e.  j Lg(j<c)  must  be  s,-135#).  Consid.r 

<*(=10,  with  present  yL^(j'a)  "“112°,  so  135°-112°“23°  more  lag  is  allowed. 

But  this  lacf  will  be  followed  by  a  lecf,  so  allow  say  10°  for  it,  giving 

23°+10°=33°  more  lag  allowable.  This  locates  the  lacf  at  15. 4  (tan33°=.65, 

and  10/ .  65=15 . 4) ,  so  we  set  the  lacf  at  <,.=  15  (i.e.  pole  at  -15),  giving 

Lg  j  (s)  =  31 . 6  ( 1  +  .  2s) _ .  Is  sketched  in  Fig.  2. 

s(l+s) (14  |^) 

Looking  ahead  at  t>'"40,  |  (j6C))j  =  -20db,  so  soon  Lg(jj)  can  make  it*; 

asymptotic  left  turn  under  the  B.^  boundary.  Our  plan  is  to  add  two  more 

leefs,  and  finally  two  complex  pole  pairs,  in  order  to  have  an  excess  ^ 

of  poles  over  zeros  of  4.  We  try  one  lecf  at  u>*40,  giving 

Lg^(s)g^‘  ^  40^  .  is  sketched  in  Fig.  2. 

s (1+s)  (1+  f^) 

We're  ready  now  for  the  last  lecf,  in  order  to  achieve  (an  asymptotic) 
horizontal  segment  for  lLg(j*o)|  ,  before  the  final  -24dl>/ octavo  slope. 

(We  foil  ow  Bode  in  this  respect,  a  good  master  to  follow.)  Whore  should 
this  horizontal  segment  be  located?  The  bottom  of  F.j  (see  Fig.  1)  is  at 
-22.5db.  Allow  2db  margin,  3db  correction  due  to  the  list  lecf,  1  .  5Jb  : 

t  I  *  *  ■  e  f  f .  t  o  t  the  1  <■  *-  I  .(I  W ,  g  i  v  i  : .  g  .<  t <  * !  .  <  i  <1 :  -  v  2  -  -  *  ■  *  ;  .  -  f  * 

We’  11  use  a  dumping  factor  of  .  6  for  the  2  complex  pole  pairs,  so  no 
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correction  need  be  allowed  for  tliem.  Thus  the  final  break,  for  |Lq(Jv)] 

asymptotic  is  to  be  at  -29db,  wliich  [Lq^(jw)|  achieves  at  u)=60.  Hence, 

the  last  leef  is  at  to“60.  The  resulting  phase  due  to  L_.  and  the  locf  at 

04 

w=60,  is  -66°.  We  could  have  -180°  at  this  point,  but  we'll  allow  an  additional 
15°  margin  (a  matter  of  taste;  it  depends  on  the  problem  —  presence  of 
higher  order  modes,  etc).  This  means  100°  phase  lag  is  permitted,  50°  due 

to  each  complex  pole  pair  (180-66-15=100).  For  C=.6,  this  locates  them  at  100. 
Thus  Ln(s)»  31.6(l+.2s)(l+fo)(H-fe)  # 


s(i.s)  (i+  Ijhl  11^4  *^,2 

Discussion 

Lq ( j(aj)  is  sketched  in  Fig.  1.  A  well-designed,  i.e.  "economical"  ( i'  .) 
is  close  to  its  boundary  B(-j)  at  each  to.  The  vertical  line  -140°  is  the 

dominating  B(io)  for  ui<5  and  the  right  side  of  B.  (line  -135°)  is  the  boundary 

n 

effectively  for  5<i<30;  so  our  Lq(.}u>)-  is  pretty  good  in  this  respect  sine.' 
it  is  pretty  close  to  these  boundaries.  There  is  tradeoff  between  complexity 
of  Lq(s)  (number  of  its  poles  and  zeros)  and  its  final  cut-off  frequency,  now 
at  uj*100.  There  is  some  phase  to  spare  between  Lq(Ju))  and  the  boundaries,  so 
use  of  more  poles  and  zeros  in  Lq(s)  would  permit  this  cut-off  frequency  to 
be  reduced  a  bit  below  100,  but  not  by  much.  On  the  other  hand,  if  ve  want  to 
reduce  the  number  of  poles  and  zeros  of  Lq(s),  we  must  pay  the  price  in  a 
larger  cut-off  frequency.  We  could  economize  significantly,  of  course, 
by  allowing  more  phase  lag  in  the  low  frequency  range.  If  -180°  was  permitted 
at  up-1,  we  could  decrease  J  Lq  ( jo»)  |  at  a  rate  of  12db/octnve;  so  with  jl.  !  = 

25db  at  <a=l,  it  would  be  13db  at  u>=2  (instead  of  the  present  16db).  Even 
with  no  more  saving,  this  5db  difference,  would  allow  a  cut-off  frequency 
at  about  70  instead  of  100. 
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Also,  Fig.  1  reveals  (immediately,  without  any  shaping  of  I.  required) 
that  reduction  (i.e.  easing)  of  the  specifications  at  ia=l  to  about  21db 
(instead  of  =26db) ,  would  have  the  same  effect  as  the  above.  One  can  check 
how  badly  the  specifications  arc  compromised  by  such  easing.  The  design 
technique  is  thus  highly  "transparent"  in  revealing  the  trade-offs  between 
performance  tolerances,  complexity  of  the  compensation,  stability  margins, 
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DESIGN  EXAMPLE:  Uncertain  Plant  P-  k 

s(s+p) 

k:  [10, HO],  p:  [-2,2]. 

Note:  For  part  of  range  of  p,  plant  is  unstable. 

-  \  «T  I’  €  1  Cm  1  ^  ^  • 


Specifications  on  T(ja») 


(see  Fig.  A) 


.  ^  .  r.  i  T 


U 

0 

.5 

1 

2 

5 

10 

EXCESS  e^ 

MAX(D3) 

0 

1.3 

1 

-6 

-17 

-33 

OVER  ZEROS 

MIN(DB) 

0 

-2 

-7 

-19 

-35 

-65 

ALSO 

(AlTl)MAX 

0 

3.3 

8 

13 

18 

32 

L _ |<3DB, 

1  +1 

i  '  l*  | 

u£.P. 

N.5.  At  high  u,  (A|t|)MAX  MUST  BE  >  |APJmAX.  OK. 


Sten  1  Calculation  and  Construction  of  Plant  Templates  { P (ju>) )~ 

The  templates  are  shown  in  Fig.  E^.  It  suffices  to  calculate  (at  a 
fixed  u;)  several  values  of  l/jwQte+'p) ,  i.e.  at  different  p  values.  This 
gives  the  bottom  curve  of  the  template.  Then  extend  vertically  by 
18db(80/10=8;  201og8=18). 

Step  2  Use  the  procedure  described  in  the  notes  (p.17  etc)  to  find  bound-, 
on  LqOoj)  in  order  to  satisfy  the  specifications  on  |T(jw)|.  A  nominal  plant 
must  be  chosen;  p*2,  k*l  was  chosen  as  nominal  and  marked  heavily  in  Fig.  B7. 
(It  helps  to  have  the  templates  on  transparent  paper  or  plastic.)  Already 
at  to*l,  the  bound  on  Lg(j),  denoted  by  B(l)  is  determined  by  |L/(1+L) 3PR, 
rather  than  by  the  constraints  on  |T(jw)|.  And  this  is  so  for  u>>l  also, 
which  Is  not  typical  for  stable  plants  but  more  likely  f~r  plants  which  ran 
be  open-loop  unstable  for  part  of  the  parameter  range. 
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Shaping  of  Lg(j<,>) 

The  boundary  at  oj”2  dominates,  i.c.  determines  the  level  of  | L^|  at 
ui“2  to  be  =10db.  The  boundary  at  .,>=5  determines  the  phase  there  to  be 
=-85°or  so.  This  necessitates  a  lead  corner  frequency  at  son..'  This 

was  chosen  to  be  at  1  and  lags  introduced  thereafter  such  that  at  u'=5  the 
phase  >  -85°;  they  were  chosen  at  5  and  8  (see  Figs.  B^.C).  A  leef  is  then 
needed  to  be  followed  by  a  complex  pole  pair  with  r,=0.6.  The  corner  below 
the  high-frequency  bound  can  be  turned  when  jL|  is  s  -25db.  If  we  try  a  leef 
at  40,  then  phase  requirements  at  (a-20  force  the  final  cut-off  frequency  to  b 
well  beyond  the  w  value  at  which  |LQ|  — 25db.  A  leef  at  u)=25  gives  compat¬ 
ibility  of  phase  needs  at  ai*=20  and  turning  the  corner  when  |l.Q|  permits  it. 
This  gives  100  as  the  frequency  at  which  the  final  complex  pole  pair  can  bo 
inserted.  This  gives  ^(s)®  4(l+s)  (1+  y-) 

s(l+.2s)(l+  g) [1+  iqo  +  ^100^ 

*GP .  -  G  (s)  ,  giving  G(s),  —  see  Figs.  C. 

s (s+2) 

Finding  F(s) 

The  proper  G  guarantees  that&|T|  does  not  exceed  tlioso  allowed.  The 
next  stop  is  to  find  Che  range  of  |L(j(l.)/(l+L(jui)  |.  Place  the  template  cf 
Che  plant  at  w=5  (for  example)'  on  the  point  Ly(j5)  in  the  Nichols  Chart, 
i.e.  on  bdb  /- 82°  and  it  is  seen  that  |L/(I+b)|  mnx=  2.4db,  min  —  1 . / J J> . 

But  the  specifications  require  |T  |=  !  F I  -  /  ( 1+1.)  ;  t  l -33 ,  ~  1 7 )  db.  Therefore,  i 
is  required  chat :  -33. 3db<  |  F|< -19. 4dh .  In  this  wav  we  obtain  the  bound.- 


on!  F(p)|  shown  in  Fig.  D.  It  is  easy  to  find  an  F(s)  which  satisfies 
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APPENDIX  B 
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Generating  the  Transfer  Function  Matrix  Models 

This  appendix  contains  the  full  state  space  matrices  from  FPCCS1M, 
the  lateral  mode  matrices  derived  by  simplification  of  the  full  state 
space  matrices,  the  transfer  function  matrices  derived  from  the 
lateral  mode  state  space  matrices  and  the  transfer  function  matrices 
generated  from  the  P^  matrices  using: 

P~ 1  ■=  £'  (3-13) 

ll  m  U/q^J  (3-14) 

The  full  state  space  matrices  (output  of  FPCCS1M)  are; 

1.  Flight  condition  1  (slow  speed  condition)  at  0.6  mach  at  sea 

level  (Eq.  B-l).  The  trim  conditions  are  0.6  mach,  the  angle  of  attack 
is  1.03  degrees,  the  dynamic  pressure  is  553.04  lbs/sq  in  and  the 
forward  velocity  is  670.5  ft/sec  (see  Table  B-l). 

2.  Flight  condition  2  (cruise  condition)  at  0.9  mach  at  30,000 

feet  (Eq.  B-2) .  The  trim  conditions  are  0.9  mach,  the  angle  of  attack 

is  1.26  degrees,  the  dynamic  pressure  is  356.15  lbs/sq  in  and  the 

forward  velocity  is  896.0  ft/sec  (see  Table  B-2). 

3.  Flight  condition  3  (high  speed  condition)  at  2.3  mach  at 
40,000  feet  (Eq.  B-3) .  The  trim  conditions  are  mach  2.3,  angle  of 
attack  is  0.63  degrees,  the  dynamic  pressure  is  1449.8  lbs/sq  in  and  the 
forward  velocity  is  2229.0  ft/sec  (see  Table  B-3). 

Simplify  the  lateral  mode  matrix  (see  Section  1V-3  for  simplifying 
assumptions  and  t:.c  matrices  become: 
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Flight  condition  1  (slow  speed  condition) 
0.6  mach  at  sea  level 
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Ax  +  Bu 
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-666.0 

13 
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-2 

'o.o 
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0. 198 
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Flight  condition  2  (cruise  condition) 


0.9  mach  at  30,000  feet 
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x  =  Ax  +  Bu 
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Flight  condition  3  (high  speed  condition) 
2.3  mach  at  40,000  feet 


x  =  Ax  +  Bu 
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Next  input  the  simplified  lateral  mode  state  space  matrices  into  a 
computer  program  to  generate  the  transfer  function  matrices,  P.  The 
resulting  nine  transfer  functions  are  of  the  form: 


_ 1 _ 

characteristics  equation  (C.E.)  of  P 


nine  numerators;  in  the  form 
ot  a  tiuco-hv-three  matrix 


8. lOOt  13  s**3  6. 253E-13  s**3 
0.0256  s**2  -0.0179  s**2 
0.0587  s**l  0.0116  s**l 
0.306  s**0  0.0692  s**0 


3. 766E-13  s**3  0.0570  s**2  6.040E-13  s**3 
0.269  s**2  -0.123  s**l  -0.0355  s**2 
0.0762  s**l  -13.568  s**0  -0.00995  s**l 
3.16  s**0  5.583  s**0 
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-35.9018  s**l  0.5527  s**l  763.5401 

0.7322  s**0  -0.01162  s**0  -15.4153 
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Eq.  B-4,  B-5  and  B-b  are  the  resulting  matrices  for  FC  III  (slow  speed 
condition),  FC  li 2  (cruise  condition)  and  FC  II 3  (high  speed  condition), 
respectively. 

Using  Eq.  3-13  and  3-14,  the  £  transfer  functions  matrices  are 
generated.  With  £  in  the  form  described  above,  the  inverse  of  £,  is  in 


the  form: 


_1  C.E.  of  P[ 3  x  3] 


determinant of  P  (de t  P) 


From  Eq.  3-14,  £  is  equal  to  the  inverse  elements  of  £  .  All  nine 
transfer  functions  have  a  set  of  common  numerators  (zeros),  the  det  £. 
Also,  all  nine  transfer  functions  have  a  common  set  of  denominators 
(poles),  the  C.E.  of  JP.  Therefore  all  nine  transfer  functions  have 
these  common  poles  and  other  poles.  The  common  poles  are  removed  from 
the  nine  terms  and  listed  separately.  Therefore  the  £  matrices  are  in 
the  form: 


common  numerator  (poles  generated  from  det  P)  3  x  3  q . 

C.E.  of  P  matrix1^ 


The  £  transfer  function  matrices  are  Eq.  B-7  (FC  It  1),  B-8  (FC  II 2)  and 
B-9  (FC  113). 

Notice  that  many  of  the  terms  have  a  singular  root  which  is  quite 
large  in  magnitude,  i.e.  >  10^  larger  than  any  other  roots.  This  very- 
far  off  root  has  very  little  effect  on  the  system  performance  until  long 
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after  the  range  of  frequencies  of  interest,  thus  it  effectively  has  no 
effect  on  the  design  except  as  a  gain.  Therefore  this  root  is  treated 
as  a  gain.  Also,  delete  the  common  roots  in  the  numerators  and  the 
denominators  (these  common  roots  are  listed  last  for  each  term  above  and 
the  simplified  matrices  are  in  Eq.  B-10  (FC  #1),  B-ll  (FC  if 2)  and  B-12 
(FC  if3). 

Simplify  and  rewrite  into  a  more  useful  form.  Thus  the  final  Q 
matrices  are: 

Flight  Condition  1  (slow  speed  condition)  0.6  mach  at  sea  level 


1 

1 

1 

0. 1747(s+0. 08111) 
(s+2.1268) 

-0.00094 (s+1 .7127) 

0. 1423(s-5. 2379) 

1 

1 

1 

=  0.03565 

0. 004 224 (s+3. 2787) 
(s-17.352) 

0.00749(s+0. 9848) 

-0.115Ks-42.2798) 

1 

1 

1 

-0.01294(s+2.0513) 

(s+46.4722) 

0. 0384 3 (s+0. 5093) 

0. 3526 (s+7 3. 6488) 

Flight  Condition  2  (cruise  condition)  0.9  mach  at  30,000  feet 
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1 


1 


1 


0.04003(s-0. 3329)  0.000499 (s+0. 6835) 

(s+1.2021) 


0.02069(s+19. 7814) 


0.005761 

0.002100(s-32.4517) 

(s+1.7274) 

0.003667 (s+0. 4 731) 

-0.05518(s-58. 9133) 

1 

1 

1 

-0.003195(s+60.919) 

(s+1.5098) 

0 . 009 1 5  5 ( s+0 .2782) 

0.08397(s+98.0476) 

Flight  condition  3  (high  speed 

condition)  2.3  mach 

at  40,000  feet 

1 

1 

1 

0. 1651(s+l .0054) 
(s-0.3343) 

0. 00 1675 (s+5 .0656) 

0.08797 (s+41 .2222) 

1 

1 

1 

0.04526 

0.007331 (s-45. 244) 
(s+1.3769) 

0.01 422 (s+0. 2047) 

-0.2126(8-149.13) 

1 

1 

1 

-0. 01 123(s+80. 350) 
(s+1.2633) 

0.03552 (s+0. 3318) 

0. 3257 (s+243. 218) 
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Appendix  C 


Calculations  and  Tables  for  the  and 

This  appendix  contains  the  calculations  necessary  to  design  and 
L^.  Use  (Table  C-2)  to  generate  the  plant  templates  for  (Fig. 
C-3).  Using: 


*33' 


-1  b22  |q„|  +  b  1 2  [  q 


1  +  L3I  < 


32' 


Ill 

l31 1 


32 


(5-2) 


the  bounds  for  each  flight  condition  for  (1  +  L^)  are  generated  (Tables 
C-2,  C-3,  C-4  and  C-7).  The  bounds  for  the  nominal  design  bounds  are 
listed  in  Table  C-8. 

Using  the  plant  templates,  transfer  the  nominal  bounds  to  the 
Nichols  Chart  to  generate  the  nominal  design  bounds  (Fig.  C-5) .  With 
the  nominal  design  bounds  drawn,  shape  the  (Fig.  C-6) .  The 

synthesized  is  in  Table  C-10.  The  compensator  for  the  third  loop,  G^> 
derived  from  L^q  is  in  Table  C— 11.  for  the  three  FC  are  listed  in 

Tables  C-10  (since  (L^) ^  was  the  nominal  case,  L^q  =  (L^)^,  C-12  and 


(L). 


_ 6.7879 _ 

s  (s  +  98 . 048 )  (s  +  0.700  +  j0.714)(s  +  0.700 


j  0 . 7 14  ■ 


(l)3 


_ 13.748 _ 

s  (s  +  243. 218) (s  +  0.700  +  j0.714)(s  +  0. 700-j0. 714) 
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Table  C- 

1  q ^  Values 

;  (magnitudes  in  decibels) 

Freq. 

FC 

in 

FC 

if2 

FC 

H2 

(rad/sec) 

(db) 

(degrees) 

(db) 

(degrees) 

(db) 

(degrees) 

0.1 

-2.5761 

-53.6500 

-9.2205 

11.9692 

-2.1774 

10.9735 

0.2 

-7.0914 

-73.2996 

-10.2728 

21.5582 

-3.2598 

19.6399 

0.4 

-12.7385 

-89.1902 

-13.1533 

31.8346 

-6.2597 

28.4176 

0.8 

-19.0523 

-104.8253 

-18.7188 

33.7688 

-12.1708 

28.6118 

1.0 

-21.2673 

-110.5459 

-21.1117 

31.8363 

-14.7268 

26.6695 

2.0 

-29.1507 

-130.9180 

-30.2719 

21.5605 

-24.3723 

17.1994 

4.0 

-38.9823 

-150.8389 

-41.2605 

11.9708 

-35.6117 

9.3316 

8.0 

-50.2375 

-164.5315 

-53.0007 

6.1633 

-47.4322 

4.7702 

10.0 

-54.0093 

-167.5286 

-56.8398 

4.9485 

-51.2815 

3.8265 

20.0 

-65.9070 

-173.6976 

-68.8307 

2.4863 

-63.2864 

1.9202 

40.0 

-77.9115 

-176.8403 

-80.8593 

1.2447 

-75.3184 

0.9610 

80.0 

-89.9435 

-178.4191 

-92.8973 

0.6225 

-87.3523 

0.4806 

100.0 

-93.8188 

-178.9352 

-96.7733 

0.4980 

-91.2335 

0.3845 

Table  C-2  q^  Values  (magnitudes  in  decibels) 


Freq. 

FC 

n 

FC 

i n 

FC 

#3 

(rad/sec) 

(db) 

(degrees) 

(db) 

(degrees) 

(db) 

(degrees) 

0.1 

-57.2600 

-0.0778 

-63.0458 

-0.0584 

-64.8543 

-0.0236 

0.2 

-57.2600 

-0. 1556 

-63.0458 

-0.1169 

-64.8543 

-0.0471 

0.4 

-57.2601 

-0.3112 

-63.0458 

-0.2337 

-64.8543 

-0.0942 

0.8 

-57.2605 

-0.6223 

-63.0461 

-0.4675 

-64 . 8544 

-0. 1685 

1.0 

-57.2608 

-0.7779 

-63.0462 

-0.5843 

-64.8544 

-0.2356 

2.0 

-57.2632 

-1.5555 

-63.0476 

-1.1686 

-64 . 8546 

-0.4711 

4.0 

-57.2728 

-3.1088 

-63.0530 

-2.3362 

-64.8555 

-0.9422 

8.0 

-57.3109 

-6. 1994 

-63.0746 

-4 . 6646 

-64.8590 

-1.8839 

10.0 

-57.3393 

-7.7323 

-63.0907 

-5.8235 

-64.8616 

-2.3544 

20.0 

-57.5690 

-15.1928 

-63.2228 

-11.5292 

-64.8836 

-4.7009 

40.0 

-58.3626 

-28.5072 

-63.7144 

-22.1937 

-64 . 9702 

-9.3393 

80. 0 

-60.6443 

-47.3670 

-65.2618 

-39.2120 

-65.3005 

-18.2072 

100.0 

-61.7987 

-53.6268 

-60 . 1 425 

-45.5648 

-65.5326 

-22.3502 

V  S  v^l 


Table  C-3 

q32  Values 

(magnitudes  in 

decibels) 

Freq. 

FC  1H 

fc  n 

FC  113 

rad/ sec) 

(db) 

(db) 

(db) 

0.1 

5.04 

6.56 

11.31 

0.2 

4.59 

5.28 

10.34 

0.4 

3.12 

2.22 

7.79 

0.8 

-0. 19 

-2.58 

3.35 

1.0 

-1.65 

-4.35 

1.65 

2.0 

-6.95 

-10.13 

-4.03 

4.0 

-12.76 

-16.09 

-9.97 

8.0 

-18.73 

-22.09 

-15.96 

10.0 

-20.66 

-24.03 

-17.90 

20.0 

-26.68 

-30.05 

-23.92 

40.0 

-32.69 

-36.06 

-29.94 

80.0 

-38.71 

-42.09 

-35.96 

100.0 

-40.65 

-44.02 

-37.89 

■£-•>* 
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Table  C-4 


q^ ^  Values  (magnitudes  in  decibels) 


•>  .V  -3 


Freq. 

fc  in 

fc  n 

FC  113 

ad/sec) 

(db) 

(db) 

(db) 

0.1 

-30.79 

-34.17 

-28.05 

0.2 

-30.82 

-34.23 

-28.13 

0.4 

-30.95 

-34.45 

-28.44 

0.8 

-31.40 

-35.23 

-29.49 

1.0 

-31.71 

-35.73 

-30.14 

2.0 

-33.69 

-38.56 

-33.47 

4.0 

-37.63 

-43.21 

-38.46 

8.0 

-43.01 

-48.86 

-44.20 

10.0 

-44.92 

-50.79 

-46.13 

20.0 

-51.35 

-57.06 

-52.29 

40.0 

-59.00 

-64.18 

-59.00 

80.0 

-68.59 

-72.99 

-67.05 

100.0 

-72.05 

-76.25 

-70.06 
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Table  05 

q^2  Values 

(magnitudes 

in  decibels) 

Freq. 

fc  m 

FC  #2 

FC  #3 

(rad/sec) 

(db) 

(db) 

(db) 

0.1 

26.89 

25.38 

14.54 

0.2 

26.85 

25.il 

14.54 

0.4 

26.68 

24.19 

14.51 

0.8 

26.05 

21.72 

14.43 

1.0 

25.65 

20.50 

14.38 

2.0 

23.17 

15.66 

13.91 

4.0 

18.81 

10.00 

12.44 

8.0 

13.32 

4.07 

9.11 

10.0 

11.45 

2.14 

7.64 

20.0 

5.53 

-3.86 

2.34 

40.0 

-0.47 

-9.88 

-3.48 

80.0 

-6.49 

-15.90 

-9.44 

100.0 

-8.42 

-17.84 

-11.38 

Table  C-6  Values  (magnitudes  in  decibels) 


Freq. 

FC  #1 

fc  n 

FC  #3 

ad/ sec) 

(db) 

(db) 

(db) 

0.1 

-26.41 

-37.03 

-38.07 

0.2 

-26.41 

-37.03 

-38.07 

0.4 

-26.43 

-37.03 

-38.07 

0.8 

-26.50 

-37.04 

-38.08 

1.0 

-26.56 

-37.04 

-38.08 

2.0 

-27.00 

-37.07 

-38.08 

4.0 

-28.40 

-37.20 

-38.12 

8.0 

-31.63 

-37.69 

-36.24 

10.0 

-33.07 

-38.02 

-38.32 

20.0 

-38.33 

-40.09 

-38.99 

40.0 

-44.14 

-44.10 

-40.96 

80.0 

-50. 10 

-49.42 

-44 . 86 

100.0 

-52.03 

-51.27 

-46.45 

07 
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Table  C-7  Boundary  magnitudes  (decibels)  for  frequencies  of  interest 


-V 


Freq. 

b22 

a22 

b32 

b  12 

So 

(rad/sec) 

(db) 

(.db) 

(db) 

(db) 

0.1 

0.04 

0.04 

-38.26 

-33.52 

% , 

0.2 

0.13 

-0.14 

-33.28 

-29.32 

0.4 

0.62 

-0.55 

-29.76 

-26.11 

r*  * 

0.8 

2.4 

-2.0 

-27.39 

-24. 14 

,\ 

1.0 

3.3 

-2.9 

-26.69 

-23.83 

2.0 

-1.1 

-8.2 

-24.98 

-23.72 

\;I 

4.0 

-9.8 

-17.1 

-25.48 

-25.06 

8.0 

-15.7 

-28.7 

-28.88 

-28.38 

10.0 

-17.3 

-32.8 

-30.42 

-29.85 

20.0 

-22.8 

-46.9 

-35.86 

-35. 16 

r- 

40.0 

-28.6 

-63.0 

-41.72 

-40.98 

80.0 

-34.5 

-80.4 

-47.70 

-46.95 

100.0 

-36.5 

-86.2 

-49.64 

-48.88 

*x 
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Table 

C-8  Loop 

bounds 

,v 

>■ 

Composite 

.V 

a; 

Freq. 

FC  in 

FC  H2 

FC  i/3 

Bound 

-S 

(rad/sec) 

(db) 

(db) 

(db) 

(db) 

S' 

0. 1 

-23.12 

-24.91 

-32.17 

-23.12 

t  • 

0.2 

-22.96 

-24.73 

-31.98 

-22.96 

1% 

0.4 

-21.76 

-23.38 

-30.37 

-21.76 

r.* 

0.8 

-17.80 

-19.00 

-25.22 

-17.80 

£ 

1.0 

-15.97 

-16.95 

-22.83 

-15.97 

2.0 

-24.51 

-25.40 

-32.38 

-24.51 

4.0 

-36.39 

-37.55 

-47.03 

-36.39 

**  , 

8.0 

-34.96 

-36.89 

-46.82 

-34.96 

*v 

10.0 

-33.09 

-34.42 

-45.18 

-33.09 

20.0 

-25.70 

-27.19 

-38.36 

-25.70 

y 

40.0 

-16.97 

-19.39 

-29.86 

-16.97 

80.0 

-8.04 

-9.25 

-19.82 

-8.04 

100.0 

-5.42 

-6.33 

-15.85 

-5.42 

C-9 


Table  C-10  L^q  Values 


10 

s(s  +  73.65)(s 

+  0.70  +  j  . 7 14) (s 

+  0.70  -  j  .714) 

Freq. 

(rad/sec) 

(db) 

(degrees) 

0.1000 

2.659b 

-98.1285 

0.2000 

-3.3624 

-106.4191 

0.4000 

-9.4659 

-124.0077 

0.8000 

-16.8167 

-162.8129 

1.0000 

-20.2668 

-180.7863 

2.0000 

-35.6310 

-228.5324 

4.0000 

-53.4861 

-252.6367 

8.0000 

-71.5782 

-266.1187 

10.0000 

-77.4215 

-269.6832 

20.0000 

-95.7139 

-281.1784 

40.0000 

-114.5895 

-296.5010 

80.0000 

-134.9131 

-316.3638 

100.0000 

-141.8820 

-322.8262 

Table  C-ll  Values 

98.9173 

s(s 

+  0.70  +  j  .714) (s  +  0.7 

-  j  .714) 

Freq. 

(rad/ sec) 

(db) 

(degrees) 

0.1000 

59.9085 

-98.0507 

0.2000 

53.8865 

-106.2635 

0.4000 

47.7831 

-123.6965 

0.8000 

40.4327 

-162.1906 

1.0000 

36.9829 

-180.0084 

2.0000 

21.6211 

-226.9769 

4.0000 

3.7756 

-249.5280 

8.0000 

-14.2783 

-259.9194 

10.0000 

-20.0933 

-261.9510 

20.0000 

-38.1560 

-265.9858 

40.0000 

-56.2181 

-267.9942 

80.0000 

-74.2799 

-268.9973 

Vs 


Appendix  D 
Calculations  for  Q„ 


This  section  contains  the  calculations  for  Q  (Chapter  VI).  For 
loop  2,  the  exact  loop  equation  is  used  to  synthesize  The  terms 

enclosed  are  used  to  generate  the  plant  template,  Q22e.  exact 

equation  is: 


n  2- 

q22[a  +  1- 

.j)(l  +  L 

3)  “  Y13^ 

Q22e 

[(1  ♦  Lx)(l 

+  L3*  "  Y13 

]  -  [y12 

a  *  l3)  *  Y23<1  *  Lj) 

-  T]  (6  3 

The  Q22e  values 

are  listed 

in  Table 

6-1. 

Table  6 

-1  Q22e 

Values 

Freq. 

FC 

//I 

FC 

#2 

FC 

in 

(rad/sec) 

(db) 

(degrees) 

(db) 

(degrees) 

(db) 

(degrees) 

0.1 

14.65 

11.98 

12.72 

21.36 

24.48 

87.27 

0.2 

15.90 

17.09 

14.62 

19.95 

26.77 

118.18 

0.4 

17.73 

18.48 

16.02 

17.45 

27.32 

-220.24 

0.8 

20.49 

12.92 

18.58 

15.97 

25.96 

-203.39 

1.0 

21.37 

5.94 

19.99 

11.33 

25.51 

-196.86 

2.0 

22.93 

-17.96 

25.61 

-23.91 

24.33 

-183.65 

4.0 

24.93 

-110.00 

15.63 

-157.93 

19.94 

-180.43 

8.0 

7.72 

-158.38 

0.076 

-168.43 

11.39 

-178.00 

10.0 

3.27 

-160.53 

-4.13 

-168.00 

8.08 

-177.15 

20.0 

-9.27 

-159.08 

-16.48 

-165.83 

-3.05 

-174.62 

40.0 

-20.64 

-148.40 

-28.13 

-156.39 

-14.75 

-170.14 

80.0 

-30.42 

-130.91 

-38.63 

-139.97 

-26.38 

-161.38 

100.0 

-33.14 

-124.94 

-41.62 

-133.75 

-30.02 

-157.26 

For  example, 

at 

FC  H  for  w 

*=  1  rad/sec;  q^  = 

10.61  db  and  -45.44 

degrees , 

a  ♦  Vj 

m 

15.88  db  and  -72.46 

degrees,  (,1 

+  Vi  -  - 

■0.89  db 

and  0.073 

degrees , 

Y 

12  “  -19-40 

db  and  - 

]  1 2 . U4  degree 

s,  Yn  =  - 

■20.23  db 

-f.  V.  S'  .-  .-.VV  V.  ; V  . 


-2.28  db  and  15.44  degrees  and  T  =  -14.74  db 


and  -94.91  degrees,  = 
and  263.3  degrees.  Using  Eq.  6-3,  solve  for  Q22e: 

0  _  10.6lZ=45.44f  (15. 88^72. 46)  (-0. 8SZa. 073)-(-20.  23/^94. 91)1 

22e  [(15.88/272.46)  (-0 . 89/0_. 0 73) - C~20 . 23^?4 . 9 1 )  ] 

-  [  (-19.40/2,112.04)  (-0.89/0. 073) +  (-2.28 ^5.44)  (15.88^72.16) 

-  (-14.74/2^3.3)] 

_  10. 6l/=45. 44 [ (14. 99Z-72.39)-(-20. 23/^94.91)1 
(14.99/^72.39)  -  (-20.23/294.91) 

-  [(-20. 29/^111. 97)  +  (13.60/257.02)  -  (-14.74^263.3)] 

.  (10.6l/r.45.44)  (15. 13/=72. 77) 

(15.13/2.72.77)  -  (4.86/257.78) 


_  25.74/^118.21 
4.16/2123.89 

^22e  ^0J=  1  rad/sec  at  FC  Ifl)  =  2 1.58/5.68 


Thus  the  nominal  bound  for  “  21.58  db  and  5.68  degrees.  After 

all  of  these  bounds  are  calculated,  then  they  are  transferred  to  the 
Nichols  Chart  to  generate  a  nominal  design  bound. 
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Appendix  E 


Designed  Controller  Responses 

This  appendix  contains  the  results  of  the  computer  simulation. 
Note  the  magnitude  scales  are  different  for  each  plot. 

The  exact  equation  was  used  to  generate  these  plots  (Eq.  7-2). 
Using  simple  algebra,  Eq.  7-2(a)  and  7-2(c)  are  solved  in  terms  of  Eq. 
7-2 (b).  The  results  are; 
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Fig.  E-6  t 


22  Responses  r(t)  =  22.4  ft/sec 
E-7 
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